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ABSTRACT 

This  thesis  provides  the  student  of  combat  modelling 
with  a  computer  program  for  a  relatively  simple  combat 
model  that  can  be  used  in  a  classroom  environment  for  study 
and  analysis.   The  model  is  an  aggregated,  force-on- force 
ground-combat  model  that  uses  Lanchester's  aimed-fire 
equations  for  casualty  assessment.   The  original  version 
of  the  model  was  developed  in  1979  in  a  previous  thesis, 
and  shortfalls  in  the  original  model  have  been  overcome, 
along  with  the  addition  of  several  enrichments  to  provide 
added  user  flexibility,   A  user's  manual  is  provided  to 
facilitate  user  access  to  the  model  from  a  permanent  disk 
in  the  V/,  R,  Church  Comuuter  Center. 
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I.   INTRODUCTION 

As  the  art  of  combat  modelling  becomes  more  advanced, 
combat  modelers  are  continually  building  more  and  more 
complicated  models.   To  the  beginning  modeler,  the  ability 
to  understand  how  these  models  operate  is  difficult,  if 
not  impossible.   As  a  student  of  combat  modelling,  I  have 
sensed  a  need  for  a  relatively  simple  model  that  could  be 
easily  studied  in  a  classroom  environment.   It  was  with 
this  in  mind  thst  this  project  was  initiated. 

At  present,  there  seems  to  be  no  simple  combat  model 
available  that  demonstrates  the  basics  of  model  building 
to  the  beginning  student.   In  1979 >  J«  Smoler  (a  student 
at  the  Naval  Postgraduate  School)  attempted  to  build  such 
a  model  for  his  thesis  research  ( Smoler, I979) ,  His  model 
was  a  deterministic,  force-on-force  computer  model  that 
used  Lanchester's  aimed- fire  equations  for  casualty  ass- 
essment. The  general  scheme  of  his  model  is  shown  in 
Figure  1. 

Although  Smoler 's  model  was  a  bold  attempt  at  a  simp- 
listic combat  model,  it  did  have  some  problem  areas  that 
warranted  investigation.   Some  of  these  problems  were 
first  discovered  during  a  class  project  in  a  combat  models 
class  (OA  ^655)  at  the  Naval  Postgraduate  School.   It  is 
the  purpose  of  this  thesis  to  study  Smoler 's  model  and  to 
undertake  an  enrichment  program  that  will  make  the  model 
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MOVEMENT 


DETECTION 


FIRE-ALLOCATION 


NO 


FIGURE  1.   GENERAL  SCHEME 


more  realistic  and  flexible,  while  maintaining  it's  trans- 
parency and  simplicity. 

Once  enhanced,  the  model  would  then  be  made  available 
for  use  as  an  instructional  tool  for  combat  modelling 
classes.   This  will  be  accomplished  through  the  develop- 
ment of  a  user's  manual  to  facilitate  the  use  of  the  model 
by  students  by  placing  the  model  on  a  permanent  disk  in 
the  W,  R,  Church  Computer  Center  where  it  will  be  easily 
accessed  by  any  desired  user. 

The  remainder  of  this  paper  will  discuss,  in  detail, 
the  problems  found  with  Smoler's  m.odel  and  the  methods 
used  to  solve  these  problems*   It  will  also  outline  some 
new  features  that  have  been  incorporated  into  the  model 
to  allow  user  flexibility.  Finally,  the  User's  Manual 
(Appendix  A)  will  provide  all  the  required  information 
to  enable  even  the  novice  user  to  utilize  the  model. 
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II.      THE  ORIGINAL  MODEL 

A.  GENERAL  DESCRIPTION 

Smoler's  original  model  of  land  combat  is  a  deter- 
ministic model  that  plays  combat  between  two  homogeneous 
forces,  a  blue  force  and  a  red  force.  The  blue  force  is 
comprised  of  three  subunits  in  a  static  defense,  with  each 
subunit  armed  with  three  TOW  antitank  missile  system.s. 
The  red  force  is  composed  of  three  subunits  of  three  tanks 
each,  attacking  on  pre-planned  routes.  The  battle  takes 
place  on  the  10  x  10  Km  Fulda  Box  that  has  been  developed 
and  used  in  the  STAR  simulation  model  (Wallace,  1978).  Since 
the  major  thrust  of  this  paper  is  to  alter  the  original 
model  to  a  form  that  will  be  easily  used  for  classroom 
instruction,  a  brief  discussion  of  the  major  components 
of  the  original  model  will  be  presented,  including  problem 
areas  that  have  motivated  model  changes. 

B.  ATTRITION  PROCESS 

The  attrition  process  in  the  original  model  utilizes 
Lanchester  "aimed-fire"  equations  with  variable  attrition 
coefficients.  The  Lanchester  equations  used  are  fairly 
simple  and  will  be  discussed  later.   However,  the  calcu- 
lation of  the  attrition  coefficients  is  of  more  immediate 
interest.   Sraoler  used  the  Eonder-Farrell  model  to  calcu- 
late the  coefficients,  A.  .,  the  rate  at  which  one  firer 
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of  unit  i  kills  unit  j  targets.  These  A^.'s  are  computed 
according  to: 

where  E(T.  .)  is  the  expected  tiae  for  one  firer  of  unit  i 
to  kill  one  target  of  unit  j.  The  E(T^.)  is  computed 
using  the  Bonder-Farrell  formula: 

E(T,  ,)  =  t   +  t.  -  t,  +  (t,  +  t,)/P(KIH)  + 
ij     a    1    n     n    I 

(t^  +  t^)/P(h|m)  X  ((l-P(hlh))/P(KlK)  + 
P(h|h)  -  P) 

where 

t  =  time  to  acquire  a  target 

t-|  =  time  to  fire  first  round  following  acquisition 
t,  =  time  to  fire  follov/ing  a  hit 
t  =  time  to  fire  following  a  miss 
t„  =  time  of  flight  of  a  round 
P  =  probability  of  a  first  round  hit 
P(h|h)  =  probability  of  a  hit  following  a  hit 
P(hlm)  =  probability  of  a  hit  following  a  miss 
P(KlH)  =  probability  of  a  kill  given  a  hit 
This  formula  holds  for  the  conditions  that  the  hit 
probability  of  any  round  depends  only  on  the  result  of  the 
previous  round  and  no  accumulated  damage  is  considered. 
Smoler  also  assumed  that  for  the  TOV/  v/eapon  system, 
P(KIH)  =  1.0  and  P(hlra)  =  P(h|h)  =  P,  thus  reducing  the 
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equation  to: 

E(T.  .)  =  t   +  t.  +  t.  +  (t^  +  t,)(l-P)/P 

Smoler  also  assumed  that  for  tanks,  P(K|H)  =  1,0  and 
t^  =  0,  reducing  the  formula  to: 

E(T,  ,)  =  t^  +  t.  +  t^(l-P)/P(hlm) 
1 J      a    -L     m 

Utilizing  these  equations  to  calculate  the  -^i-j's,  the 

attrition  during  each  time  step  was  computed  using  the 

Euler-Cauchy  differencing  equations  to  approximate  Lan- 

chester's  fo^ce-on-force  attrition  differential  equations, 

C.   UNIT  LOCATION  AND  MOVEMENT 

Red  and  blue  unit  locations  v/ere  handled  in  two  dif- 
ferent v;ays  in  the  original  model.   The  blue  locations 
were  left  as  user  inputs,  while  the  red  locations  had 
been  pre-determined  by  the  model  builder  and  could  not  be 
altered  by  the  user.  This  allowed  flexibility  of  defensive 
positions,  but  because  of  the  method  of  determining  movement 
routes  for  the  attackers,  this  flexibility  was  limited. 

The  method  utilized  for  route  determination  was,  at 
best,  unrealistic.   For  each  original  red  location,  a 
straight  west  to  east  route  was  calculated.   Each  route 
was  divided  into  40  meter  intervals,  since  each  red  unit 
was  assumed  to  move  that  distance  during  each  10  second 
time  interval.   This  approxim.ated  an  average  rate  of  move- 
ment of  9  m.p.h.  This  method  of  route  selection,  although 
easy  to  implement,  has  several  significant  shortcomings. 
It  does  not  allow  the  attacking  units  to  utilize  terrain 
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features  during  movement.   Also,  it  permits  attackers  to 
move  over  terrain  that  is,  in  real  life,  impassible, 

D.   BATTLE  TEP.MINATION 

In  any  combat  model,  adequate  battle  termination  rules 
must  be  considered,   Smoler  utilized  two  criteria  to 
terminate  the  battle.   The  first  of  these  was  annihilation 
of  one  of  the  two  forces.  This  criteria  is  reasonable  in 
an  expected  value  model  like  this  and  was  adequately 
handled.  This  is  not  the  case  for  the  second  termination 
rule. 

The  second  rule  for  termination  is  that  the  distance 
between  red  and  blue  forces  becomes  too  small.  The  use 
of  this  criteria  for  a  stopping  rule  is  fine,  provided  it 
is  implemented  properly.   In  the  original  model,  a  geographic 
center  was  calculated  for  each  force.   It  was  these  center 
of  mass  points  that  were  compared  to  determine  if  the 
units  were  too  close.  This  method  can  lead  to  problems, 
such  as  those  that  v/ere  encountered  during  some  initial 
trial  runs.   In  one  battle,  the  red  units,  during  their 
advance,  had  eliminated  one  entire  blue  unit  and  had  suf- 
fered the  loss  of  one  of  their  own  units.   Since  the  red 
units  are  allowed  to  advance  as  much  as  I50  meters  apart 
and  blue  units  were  set  up  so  they  were  approximately 
1300  meters  apart,  there  was  considerable  distance  between 
units  on  the  battlefield.   Since  the  two  units  had  been 
destroyed,  the  center  of  mass  computation  allowed  the  red 
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units  to  move  completely  through  the  blue  defense,  while 
both  sides  still  had  forces  remaining  to  fight  with.   This 
reflects  a  problem  with  the  center  of  mass  computation  as 
a  method  for  determining  distances  between  units. 

E.   DETECTION  AND  FIKE  ALLOCATION 

Both  detection  and  fire  allocation  processes  are  hand- 
led well  in  the  original  model  and  no  changes  have  been 
made,   A  detailed  description  of  these  methods  is  contained 
in  Smoler's  Thesis (Smoler, 1979) •  However,  a  brief  and 
general  discussion  of  both  processes  ^^11  be  included  here. 

The  detection  phenomena  is  modeled  in  tv;o  ways.  First, 
a  non- firing  detection  can  occur  as  a  result  of  an  observer's 
random  search  within  his  designated  section  of  responsibility, 
A  30  field  of  vision  for  an  observer  is  assumed,  and  the 
probability  that  an  observer  is  looking  in  the  direction 
which  enables  him  to  detect  a  target  is  computed  by  inte- 
grating the  Limicon  Function  over  limits  that  are  ±  15° 
from  the  primary  direction  the  observer  is  looking.   The 
Limicon  Function,  f(9),  is  the  following  probability  den- 
sity function: 

f(e)  =  A  +  B  X  cos(e) 

where 

D  =  assigned  sector  width/2 
A  =  -B  X  cos(D) 

B  =  1  /  2(sin(D)-  D  x  cos(D)) 

e  =  primary  direction  observer  is  looking 
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Also,  A  and  B  are   chosen  such  that 

J     f(e)de  =  1 

-D 

The  second  method  of  detection  is  a  firing  detection. 
This  phenomena  occurs  when  the  following  happens.  If  a 
firer's  location  is  within  ±  15  of  an  observer's  primary 
direction  of  observation  when  he  is  firing,  he  is  assumed 
to  be  detected  and  is  added  to  the  observer's  target  list. 
This  models  the  detection  of  a  firer  by  locating  a  firing 
signature  of  a  weapon. 

The  fire  allocation  process  is  also  modeled  well  by 
Sraoler.   Since  each 'firing  unit  is  not  restricted  to  firing 
at  one  target,  a  fraction  of  each  firing  unit  is  allocated 
to  each  target  on  the  firer's  target  list.   This  fraction 
is  determined  as  a  function  of  range  to  each  target  and  of 
predetermined  priorities  of  fire  that  are  set  so  each  unit's 
fire  priority  goes  to  those  targets  to  his  immediate  front. 
Targets  to  the  firer's  flanks  are  then  allocated  a  smaller 
percent  of  the  available  firepower.   Again,  the  detection 
and  fire  allocation  processes  have  been  handled  well  and 
no  changes  to  these  routines  v/ill  be  made. 
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III.      MODEL  CHANGES 

A.  GENERAL 

As  the  previous  section  outlines,  there  were  several 
major  problems  discovered  in  Smoler's  original  model.   In 
order  to  make  the  model  more  flexible  and  classroom  useable, 
several  major  changes  were  found  to  be  necessary.   This 
chapter  will  provide  a  detailed  description  of  these 
changes,  which  include  conditionally-deterministic  attrition, 
unit  locations  and  movement,  battle  termination  and  move- 
ment of  defenders. 

B.  CONDITIONALLY-DETERMINISTIC    ATTRITION 

In  order  to  introduce  a  stochastic  process  into  the 
attrition  computation  for  this  model,  several  options  were 
available.   These  include  the  use  of  a  Markov-process  to 
determine  casualties  or  the  use  of  random  attrition-rate 
coefficients.   After  considering  both  options,  and  in 
keeping  with  the  transparent  nature  of  this  model,  it  v/as 
decided  that  stochastic  attrition  coefficients  should  be 
utilized,  because  the  attrition-rate  coefficient  is  a 
random  quantity  measuring  a  unit's  fighting  ability  and 
can  be  realized  before  any  given  battle.   The  following 
procedure  was  used  to  develop  a  method  for  random  attrition- 
rate  coefficients. 

The  attrition-rate  coefficient,  A.  .,  is  a  measure  of 
the  rate  a  firer  of  unit  i  attrits  a  target  of  unit  j. 
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This  also  can  be  interpreted  as  a  measure  of  the  fighting 
ability  of  an  element  of  type  i.   It  is  intuitively  obvious 
that  this  is  a  variable  quantity  that  can  be  affected  by 
many  different  factors,  including  weather,  espirit-de-corps, 
previous  engagements  and  leadership,  to  mention  only  a  few. 
In  a  homogeneous  model,  one  can  attempt  to  capture  these 
random  effects  by  developing  a  distribution  of  initial 
fighting  unit  capabilities,  i.e.,  initial  A.  .'s  for  each 
unit.   Since  no  data  is  readily  available  that  captures  this 
phenomena,  an  attempt  was  made  to  fit  a  distribution  to 
an  intuitive  feeling  as  to  how  this  fighting  ability  varies 
from  one  unit  to  the  next  prior  to  a  given  engagement.   It 
was  the  author's  intuitive  feeling  that  the  A .  . ' s  should  be 
distributed  between  ,3  and  .8  with  the  majority  of  the  units 
being  rated  between  ,5  and  .6. 

Utilizing  the  above  intuition,  an  attempt  v/as  m.ade  to 
"fit"  a  distribution  to  these  assumed  fighting  levels. 
Graphically,  the  distribution  would  look  like  Figure  2. 


Figure  2.   Distribution  of  Initial  A-  .'s 

1 . 1 
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Initially,  an  attempt  was  made  to  fit  a  Beta  distribution 
to  this  curve.   The  fit  was  fairly  accurate  in  the  middle 
range,  but  was  not  satisfactory  in  the  tails.   The  same 
was  true  for  a  Normal  approximation.   It  was  finally 
determined  that  the  best  fit  was  a  straight  quadratic  fit, 
v/here  a  Uniform(0,l)  input  variable  could  be  used  and  the 
output  would  yield  the  desired  random  number.   The  quad- 
ratic fit  that  was  used  is: 

A°  =   -2U^  +  2U  +  .3 

Once  a  distribution  for  the  attrition-rate  coefficients 
had  been  derived,  the  implementation  into  the  model  was 
accomplished.   Since  it  was  assumed  that  the  fighting 
ability  of  each  unit  is  a  random  quantity  prior  to  a  given 
battle,  it  v/as  only  necessary  to  obtain  a  realization  of 
the  random  variable  for  each  unit  prior  to  the  initiation 
of  the  battle.   This  realization,  A.  .,  is  determined  by 
using  a  random  Uniforra(0,l)  number  and  the  above  formula. 
A  new  attrition-rate  coefficient  for  each  unit  is  then 
computed  during  each  time  step  using  the  equation: 


?/here 

r   =  maximum  effective  range  of  a  firer's  weapon 
r  =  current  range  between  firer  and  target 
A.  .  =  realization  of  unit's  fighting  ability 
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This  equation  is  utilized  because  it  varies  with  range, 
but  it  also  is  a  function  of  aV  ,  thus  creating  a  different 
attrition-rate  curve  for  each  unit,  depending  on  that  unit's 
fighting  ability  prior  to  the  battle.   Graphically,  as  an 


example,  the  attrition-rate  coefficient  curve  for  an  A^^  . 
of  ,6  and  a  maximum  effect: 
would  look  as  in  Figure  3, 


of  ,6  and  a  maximum  effective  range,  r^,  of  3000  meters 


A 


ij 


1000 


2000    3000 


range(meters) 


Figure  3.   Attrition-rate  coefficient  curve  for 
A?.  =  ,6  and  r^  =  3000  ra. 

Once  the  above  method  for  determining  attrition-rate 
coefficients  had  been  selected,  it  was  coded  and  included 
as  a  user  option  in  the  program.   Once  the  code  had  been 
implemented,  several  runs  were  made  utilizing  different 
random  number  seeds  to  compare  battle  results  using  this 
method  for  attrition-rate  coefficient  determination.   These 
runs  were  then  compared  with  a  run  using  the  original 
deterministic  method.   The  results  are  summarized  below. 

First,  the  model  was  run  using  the  deterministic  attrit- 
ion method.   For  purposes  of  comparison,  all  locations. 
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rnoveraent  rates,  force  levels  and  other  input  variables  were 

held  constant  for  all  runs.   Combatant  attrition  using  the 

deterministic  process  as  a  function  of  time  is  shown  in 
Figure  k* 


Force 
Levels 


Red 

Blue 

■  -  -  -  2.84 
1.01 


00   200   300   400   500 
Time 


Figure  4,   Attrition  using  deterministic  attrition 

In  this  battle,  termination  occurred  as  a  result  of  forces 
being  too  close  together,  with  force  levels  at  termination 
reflected  in  Figure  Zf   This  was  then  compared  to  the  runs 
that  utilized  the  new  attrition  module.   The  results  of 
four  of  these  runs  are  in  Figure  5,   The  A?  . ' s  are  the 
realizations  of  the  random  variable  denoting  a  unit's 
initial  fighting  capability  prior  to  the  battle.  The  three 
numbers  under  the  red  and  blue  headings  are  the  realizations 
for  each  subunit  in  the  battle  on  which  casuality  assessments 
are  computed. 

These  results  show,  in  fact,  that  a  randomly  selected 
A.  .  does  have  an  effect  on  the  final  outcome  of  the  battle. 
This  module  now  provides  the  user  the  option  of  selecting 
a  conditionally-deterministic  process  that  still  uses 
Lanchester's  equations  as  the  basis  for  the  attrition 
computation. 
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A° 

Force 

Level 

Blue 

Red 

Blue 

Red 

Time 

Result 

.33  .77  .46 

.79  .49 

.56 

1.7 

1.5 

780 

Too  close 

.39  .67  .58 

.63  .78 

.36 

4.3 

0 

330 

Blue  vrlns 

.78  .67  .80 

.74  .68 

.72 

1)*^ 

0 

320 

Blue  wins 

.32  .43  .68 

.51  .77 

.65 

1.2 

3.8 

640 

Too  close 

Deterministic 

1.1 

2.8 

460 

Too  close 

Figure  5.   Comparison  of  Stochastic  and  Deterministic 
Battles 


C.   UNIT  LOCATION  AiND  MOVEMENT 

As  was  pointed  out  in  an  earlier  section,  the  unit 
locations  and  movement  logic  that  Smoler  used  tended  to 
be  rigid  and  inflexible.   The  changes  that  have  been  made 
to  rectify  this  and  allow  more  flexibility  for  the  user 
will  be  discussed  here. 

First,  the  user  is  now  responsible  for  inputing  all 
combatant  locations,  both  attackers  and  defenders.  The 
format  for  these  inputs  will  be  specified  in  Appendix  A, 
the  User's  Manual.   Aditionally,  the  user  is  given  the 
option  of  selecting  routes  for  the  red  advance  or  allowing 
the  model  to  calculate  straight  west  to  east  routes  as  in 
the  past.  To  add  more  flexibility,  the  user  v/ill  also  be 
required  to  input  the  rate  of  advance  (vehicle  speed)  for 
the  attacking  forces.   The  inputed  rate  of  advance  is 
used  in  both  methods  of  route  calculation,  thus  controlling 
the~*red  attack. 
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1.  Straig:ht  Line  Routes 

The  option  to  use  straight  line  routes  from  west 
to  east  is  handled  in  the  same  manner  as  before.  The  route 
is  broken  down  into  discrete  distance  intervals  from  the 
initial  red  location  straight  to  the  east.  The  only  dif- 
ference is  that  the  length  of  each  interval  is  determined 
by  the  vehicle  speed  that  the  user  inputs.  The  route  for 
each  red  unit  consists  of  125  equal  length  intervals. 

2.  User  Determined  Routes 

The  option  to  allow  the  user  to  select  routes  for 
the  red  advance  has  been  added  to  enable  the  user  to  make 
use  of  available  terrain  and  to  add  realism  to  the  model. 
The  method  for  calculating  these  routes  is  straightforward 
and  is  discussed  next. 

The  user  is  required  to  input  the  original  location 
of  each  red  subunit  and  from  one  to  ten  nodes  he  wishes 
each  attacking  subunit  to  move  through.   This  information, 
along  with  vehicle  speed,  is  used  to  calculate  route  inter- 
vals that  move  the  attacking  unit  through  each  of  the 
designated  nodes.   A  complete  route  would  look  like  the  one 
depicted  in  Figure  6.  The  method  used  to  compute  the  routes 
is  as  follows. 

The  straight  line  ground  distance  between  the  first 
two  adjacent  nodes,  DIST,  is  calculated  as  shown  in  Figure 
?•   The  angle  between  the  desired  direction  of  movement 
and  straight  west  to  east  movement,  a,  is  then  calculated. 
Utilizing  these  quantities  and  the  distance  desired  to 

21 


03 
0) 
4-> 

o 
a 

c 

•H 

e 

CD 
-P 
0) 

Q 

t-i 
0) 
U3 

1=3 


•H 


22 


move  during  each  time  step,  DST,  the  distance  to  be  moved 
in  the  x  and  y  direction,  XLIT  and  YLII,  is  now  computed  as 
shown  in  Figure  7,   These  distances  are  then  added  to  the 
coordinates  of  the  previous  interval  endpoint,  point  C 
in  the  figure,  to  determine  the  coordinates  of  the  next 
interval  endpoint,  point  D,   This  same  distance  is  again 
added  to  compute  the  coordinates  of  the  next  endpoint, 
point  E,   This  process  is  continued  until  the  distance 
from  the  last  endpoint  computed  to  the  next  node  is  less 
than  DST.   This  process  is  repeated  between  the  next  tv/o 
nodes  until  the  entire  route  is  completed.   In  order  to 
insure  all  intervals  are  of  equal  length,  the  computation 
of  the  first  interval  between  any  two  nodes  must  be  con- 
sidered separately,  by  taking  into  account  the  distance 
left  over  from  the  last  computation  between  the  previous 
two  nodes.   To  do  this,  the  first  interval  takes  the  remi- 
aining  distance,  e,  and  adds  it  to  an  interval  length  of 
DST-e  for  the  first  interval  between  any  tv;o  nodes.   This 
insures  that  each  interval  along  the  route  is  of  length  DST, 
which  is  the  required  length. 

D.   BATTLE  TERMINATION 

As  previously  stated,  the  battle  is  terminated  either 
by  annihilation  of  one  of  the  two  forces  or  by  forces 
getting  too  close.   The  annihilation  criteria  has  not  been 
changed.   However,  to  insure  the  distance  criteria  is 
effective,  two  changes  have  been  made. 

First,  the  distance  between  units  is  no  longer  calculated 
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NODE        COORDINATES 

A  XIC(N-1 ),YIC(N-1 ) 

B  XLOC(I,J) 

C  XIC(N),YIC(N) 

D  XIC(N+1 ),YIC(N+1 ) 

E  XIC(N+2),YIC(N+2) 

F  XL0C(I,J+1) 


USE 

(N-1)^^  Interval  Endpoint 

User  Inputed  Node 

N    Interval  Endpoint 
/.T^^xSt  Interval  Endpoint 

/j^^pxnd  Interval  Endpoint 
User  Inputed  Node 


DIST  =   NXL^  +  YL^ 

a  =  Tan"^  (Y/XL),  where  Y  =  |yl| 

e  =  distance  less  than  DST  at  end  of  calculation 

of  intervals  between  adjacent  nodes. 
YU]   =  DST  X  SIN(a)        XLN  =  DST  x  COS(a) 


Figure  7.   Route  Computation 


^k 


as  a  center  of  mass  distance.   Instead,  a  distance  is 
computed  between  each  attacking  subunit  on  which  casualties 
are  being  computed  that  is  still  alive  and  each  defending 
subunit  that  is  still  in  the  battle.   If  any  of  these 
distances  becomes  too  small,  the  battle  is  considered  to 
have  moved  to  a  "close-in"  combat  (hand-to-hand)  that  this 
program  does  not  currently  model.   For  this  reason,  the 
battle  is  then  terminated.   However,  to  insure  the  attackers 
do  not  pass  by  the  flanks  of  remaining  defenders  and  remain 
outside  termination  distance,  a  check  of  x  coordinates 
for  each  subunit  is  also  m.ade.   If  any  attacking  subunit 's 
X  coordinate  places  him  past  the  location  of  the  forward 
most  blue  defensive  subunit  still  in  the  battle,  the  battle 
is  also  terminated. 

The  criteria  for  being  too  close  has  been  left  as  a 
user  determination.   It  is  one  of  the  inputs  that  is 
described  in  the  User's  Manual,  Appendix  A.   This  allows 
for  flexibility  of  breakpoint  distances  which  also  lends 
itself  to  the  study  of  optimum  breaking  distances  for 
various  weapon  system.s  on  the  battlefield. 

E.   ALTEI?NATE  DEFENSIVE  POSITIONS 

The  last  addition  to  the  model,  added  to  increase  the 
flexibility  and  realism  of  the  battle,  is  the  option  of 
having  the  defending  units  move  to  alternate  positions  if 
the  attackers  close  to  within  a  user  specified  range. 
The  move  is  handled  in  a  simple  and  transparent  manner. 

If  the  user  selects  this  option,  the  defender  will  move 
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v/hen  the  attacking  forces  close  within  the  breakpoint 
distance,   7/hen  this  occurs,  each  defensive  unit  that  is 
still  alive  v/ill  move  to  an  alternate  position  the  user 
has  selected.   The  duration  of  the  move  is  also  a  user 
input.   He  simply  specifies  the  number  of  10  second  time 
steps  he  wants  to  allow  the  blue  forces  to  move.   At  the 
completion  of  the  move,  during  which  the  red  forces  continue 
to  attack,  the  battle  will  continue  until  one  of  the  term- 
ination criteria  discussed  above  has  been  met. 
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IV.   FUTURE  ENHAriCSMSIIT 

Although  the  above  changes  have  been  implemented  and 
a  more  useable  and  flexible  model  has  been  created,  there 
are  several  areas  that  could  still  be  considered  in  future 
v/ork.  These  include  the  introduction  of  heterogeneous 
forces,  explicit  computation  of  ammunition  expenditure 
and  artillery  (indirect  fire).   Possible  methods  for 
employing  these  ideas  are  discussed  below. 

A.   HETEROGEriEOUS  FORCES 

- —  The  current  model  involves  combat  between  homogeneous 
forces  only.   In  other  words,  each  opposing  force  is 
comprised  entirely  of  one  weapon  system  type.   Added 
flexibility  could  be  attained  by  allowing  heterogeneous 
force  structures  on  both  sides.   This  would  enable  the  user 
to  investigate  the  effect  of  different  force  mixes  on 
battle  outcome. 

The  introduction  of  different  weapon  systems  within 
a  single  unit  would  require  extensive  changes  to  the  att- 
rition process.   Although  Lanchester  equations  for  aimed- 
fire  could  still  be  utilized,  casualty  assesment  against 
each  weapon  system  type  by  each  opposing  weapon  system 
type  v/ould  have  to  be  calculated.   The  total  attrition  of 
any  particular  unit  would  then  be  the  sum  of  the  damage 
assessed  to  each  weapon  system  of  that  unit.   Aditionally, 
the  force  level  for  each  weapon  system  would  be  required 
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as  an  input  as  well  as  a  separate  set  of  hit  and  kill 
probabilities  for  each  weapon  system  type  against  each 
type  of  target. 

Introduction  of  this  feature  would  create  a  more  comp- 
licated model,  adding  realism  but  detracting  from  the  current 
simple  and  transparent  form.   Since  the  purpose  of  the 
current  effort  has  been  to  m.aintain  this  simplicity,  this 
option  has  not  been  included,  but  could  be  considered  for 
future  work. 

B.   AiMMUNITICN  EXPENDITURE 

The  combat  process  is  a  complicated  and  intricate 
process.  One  of  the  hardest,  and  not  many  times  attempted, 
areas  to  model  is  the  logistic  area.   Hov/ever,  in  the 
model  currently  being  studied,  the  problem  of  ammunition 
expenditure  could  be  modeled. 

Since  the  model  is  an  expected  value  type  model,  it 
would  not  be  effective  to  model  the  amount  of  ammunition 
expended  by  simply  counting  bullets  as  they  are  fired. 
Actually,  in  this  type  of  model  it  would  not  be  possible 
to  count  each  round,  because  of  the  aggregated  nature  of 
each  unit.   The  only  way  to  model  ammunition  useage  is  to 
model  it  in  the  same  way  as  attrition.   This  simply  means 
that  each  unit  would  have  to  be  given  a  starting  level  of 
ammunition  on  hand,  a  basic  load,  and  each  time  the  unit 
fired,  the  expected  ammount  of  ammunition  expended  would 
be  subtracted  from  what  is  on  hand.  The  ammount  expended 
would  necessarily  be  a  function  of  the  size  of  the  unit 
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firing,  the  number  of  targets  fired  at  and  the  rate  this 
unit  is  firing.   This  can  be  directly  correlated  to  the 
attrition-rate  coefficient  for  the  firing  unit,  ':Vhen 
ammunition  on  hand  reaches  the  zero  level,  the  unit  v/ould 
have  to  be  removed  from  the  battle  permanently  or  for  some 
specified  time  period,  to  simulate  resupply. 

C.   ARTILLERY 

In  any  armed  conflict,  there  are  more  ways  to  inflict 
casualties  than  just  direct  fire.   Many  types  of  indirect 
fire  are  utilized  on  the  battlefield,  including  artillery, 
close  air  support  and  naval  gunfire.   To  include  this  type 
of  play  in  this  model,  primarily  artillery,  v;ould  involve 
a  major  effort  with  large  program 'additions. 

To  model  artillery  in  an  expected  value  model  that 
already  uses  Lanchester  equations  for  direct  fire  would 
suggest  the  use  of  Lanchester 's  area-fire  equations.  These 
equations  are  no  more  complicated  to  handle  than  the  aimed- 
fire  equations  that  are  already  in  use.   However,  other 
model  considerations  would  have  to  be  investigated.  These 
would  include  location  of  artillery  units,  v/hether  to  use 
forward  observers  for  target  location  or  only  use  pre- 
planned strikes  or  both  and  fire  allocation  procedures 
to  be  used.   Also,  a  significant  ammount  of  data  would  need 
to  be  collected  concerning  weapon  types  to  be  used,  eff- 
ective ranges,  killing  radii  and  accuracy  data. 
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To  model  artillery  would  be  a  significant  addition  to 
the  model.   However,  this  too  would  effect  the  simplicity 
and  transparency  the  model  currently  possesses*   This 
addition  has  also  been  left  for  future  consideration. 
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APPENDIX  A 
User's  Manual 

I.   INTPODUCTION 

Smoler's  model  of  land  combat  is  a  force-on- force 
combat  model  that  utilizes  Lanchester's  aimed  fire  equations 
for  casualty  assessment.  The  battle  is  simulated  on  a 
10  X  10  KM  piece  of  terrain  representing  an  area  east  of 
Fulda,  V/est  Germany,   A  portion  of  the  terrain  map  is 
inclosed  as  Figure  1.   It  is  a  computer  model  that  is 
coded  in  FORTRAN  and  is  available  for  use  in  the  '.V.  R, 
Church  Computer  Center, 

The  purpose  of  this  manual  is  to  familiarize  the  user 
with  the  model  and  to  provide  the  required  Job  Control 
Language  (JCL)  and  user  inputs  to  run  the  model.   A 
sample  listing  of  program  output  will  also  be  provided  to 
familiarize  the  user  with  the  expected  program  output. 
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II.   PPCGI?AI-I  STRUCTURE 

Smoler's  enhanced  model  is  a  computerized  model,  coded 
in  FCRTRAII,  containing  a  main  program  and  nine  subroutines. 
To  assist  the  user  in  understanding  hov;  the  model  operates, 
a  brief  description  of  the  function  of  each  of  these  major 
parts  of  the  program  will  be  included. 

A.  MAIN  PRCGRAM 

The  main  program  has  several  important  functions. 
All  of  the  input  and  output  functions,  except  the  line-of- 
sight  data,  are   contained  in  the  main  program.   Aditionally, 
the  main  program  is  used  to  structure  all  of  the  other 
functions  during  each  10  second  time  step.   Attrition, 
detection,  movement  and  fire  allocation  are  also  handled 
in  the  main  program.   The  nine  subroutines  provide  needed 
input  numbers  for  the  above  calculations.   The  general  flow 
of  the  program  is  depicted  in  the  flow  chart  shown  at 
Figure  2. 

B.  SUBROUTINES 

There  are  nine  subroutines  included  in  the  model.   Each 
of  these  subroutines  performs  a  distinct  function,  each  of 
which  will  be  discussed  below. 

1.   Subroutine  SETUP 

This  subroutine  is  used  to  read  in  the  terrain 
data  for  the  Eulda  Map.   This  terrain  data  will  be  used 
when  computing  line-of-sight  between  targets  and  observers, 
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as  well  as  providing  a  grid  system  for  unit  locations  and 
movement. 

2.  Subroutine  JXS 

This  module,  developed  by  Professor  Jim  Hartman 
of  the  Naval  Postgraduate  School,  computes  a  percent  of  a 
target  visible  to  a  particular  observer,  given  the  coor- 
dinates (location)  of  both.  This  visible  fraction  is  used 
in  the  detection  and  attrition  modules  in  the  model. 

3.  Subroutine  KCVSP 

This  subroutine  is  used  by  subroutine  IDS   in  deter- 
mining what  portion  of  a  particular  target  is  covered  by 
the  terrain  between  the  target  and  observer.   This  number 
is  used  in  the  detection  and  attrition  modules, 
if.   Subroutine  ETK 

This  module  computes  the  expected  time  for  a  given 
firer  to  kill  a  given  target.  The  calculation  is  a  function 
of  range,  time  of  flight  for  a  round  and  hit  and  kill 
probabilities  for  the  firing  weapon  system.   It  is  a  number 
that  is  used  in  the  comiputation  of  the  deterministic  attrition 
coefficients,  A.  .. 

5.   Subroutine  STOCH 

This  is  used  when  a  user  has  selected  a  stochastic 
attrition  option  to  compute  the  attrition  coefficients 
during  each  time  step.  The  calculation  is  a  function  of 
the  original  stochastically  determined  A •  .  as  well  as  a 
function  of  range. 
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7.  Subroutine    LAT-IBDA 

This  subroutine  is  used  in  conjunction  with  the 
LOS  routine  to  compute  the  detection  rate  of  a  target  by 
an  observer  given  the  percent  of  target  visible  to  the 
observer. 

8.  Subroutine  ROUTE 

This  is  used  to  compute  the  route  of  each  attacking 
red  unit  when  the  user  has  selected  the  option  of  inputing 
attacker  routes.   It  calculates  the  coordinates  of  each 
interval  endpoint  along  the  route,  making  each  interval 
length  (distance  moved  during  a  10  second  time  step)  the 
same.   The  interval  length  is  determined  by  the  speed  the 
user  has  selected  and  inputed  for  the  current  battle. 

9.  Subroutine  SLSV 

A  subroutine  that  is  used  to  calculate  the  terrain 
elevation  for  a  given  set  of  X,  Y  coordinates.   This  function 
is  used  in  conjunction  with  the  LOS  subroutine  in  computing 
line-of-sight  between  observer  and  target. 
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III.   AVAILABI-^  OPTIONS 

The  enhanced  version  of  this  model  has  been  v/ritten 
to  allow  for  maximum  user  flexibility  v/hile  maintaining 
the  simple  and  transparent  nature  of  the  model.   To  allow 
this  flexibility,  there  have  been  several  user  options 
incorporated  into  the  m.odel.   Each  of  these  options, 
including  user  responsibilities,  will  be  discussed  here, 
with  the  required  input  data  for  each  being  outlined  in 
the  next  section, 

A.   STOCHASTIC  VS  DSTEI^MIMISTIC  ATTRITION 

The  user  is  required  to  specify  whether  he  wants 
to  use  stochastic  or  deterministic  attrition  calculations. 
Both  m.ethods  utilize  Lanchester's  aim.ed-fire  equations, 
the  difference  being  the  method  of  calculating  the  attrition- 
rate  coefficients.   The  deterministic  procedure  uses  the 
Eonder-Farrell  method  of  calculating  the  A ■  . ,  while  the 
stochastic  m^ethod  selects  an  initial  random  A .  .  for  each 
unit,  and  uses  this  as  a  function  of  range  to  calculate 
the  A.  .  for  each  time  step. 

E.   ATTACK  ROUTES 

The  second  major  option  available  to  the  user  is 
the  method  of  route  comiputation  for  the  attacking  forces. 
The  user  has  to  decide  whether  to  allow  the  program  to 
compute  straight  west  to  east  routes  or  to  input  the  routes 
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he  desires  each  attacking  unit  to  follow.   To  select  his 
own  routes,  a  user  must  input  the  number  of  nodes  he 
;vishes  to  have  on  each  route  and  the  coordinates  of  each 
of  these  nodes.  The  program  will  then  compute  routes 
through  each  node.  The  nodes  must  be  inputed  in  order 
from  west  to  east  and  should  not  create  an  angle  between 
the  west  to  east  axis  and  the  route  direction  that  exceeds 

45°. 

C.  ALTERNATE   DEFEIISIVS   POSITIONS 

The  third  option  the  user  must  consider  is  the 
use  of  alternative  defensive  locations.   This  option  permits 
the  user,  if  he  desires,  the  capability  to  move  the  defen- 
ders to  alternate  positions  if  the  attackers  close  within 
some  specified  distance.  This  breakpoint  distance  is 
decided  and  inputed  by  the  user  and  is  also  used  as  the 
distance  for  battle  termination.   The  alternative  to  moving 
the  defenders  is  to  terminate  the  battle  when  the  breakpoint 
distance  is  reached. 

D.  OTHER  INPUTS 

There  are  other  inputs  that  are  required  by  the 
program  that  the  user  must  provide.  These  include  force 
sizes,  weapon  characteristics,  unit  locations  and  hit  and 
kill  probabilities.  The  required  formats  for  all  inputs 
are  outlined  in  the  next  section. 
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IV.   INPUT  DATA 

All  user  input  data  is  read  in  from  cards  at  the  start 
of  the  main  program,  MAIN,   A  brief  description  of  each 
data  element  as  v/ell  as  the  required  format  for  inputing 
this  data  is  discussed  next.   A  sample  data  deck  has 
also  been  included  as  Figure  3)  with  the  referenced  data 
element  number  appended  on  the  right.  This  sample  data 
deck  reflects  the  use  of  all  data  elements.   However, 
not  all  the  elements  described  below  are  required,  and 
those  that  are   needed  only  if  a  particular  option  is 
selected  are  noted  with  an  asteric  next  to  che  data  element 
number. 

Data  Element Description  and  Format 

1,  Selection  of  attrition  option  and  random 

number  seed,   A  1  is  entered  for  deter- 
ministic or  a  0  for  stochastic  attrition. 
This  number  is  followed  by  a  five  digit 
random  number  seed.  Format  for  this  card 
is  II, IX, 15. 

2,  Number  of  blue  units  (^IBU)  and  number 

of  red  units  (riRU),   Both  of  these 
elements  are  two  digit  integers.   Card 
format  is  I2,iX,I2. 

3»  Effective  v/eapon  ranges  listed  in  the 

following  order.  Red  minimum,  red 
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maximum,  blue  minimum  and  blue  maximum. 
Card  format  is  /f(F6.1,lX). 
Z^.,  Red  force  levels  for  each  unit. 

Card  format  is  F3.1,  IX.  The  card  v/ill 
contain  the  number  of  entries  that  equals 
the  number  of  red  units  (I^IRU)  in  data 
element  2. 

5.  Type  of  route  com.putation  desired  followed 

by  vehicle  speed.   Use  a  1  for  user 
determined  routes  or  a  0  for  program 
determined  routes.   A  one  digit  entry  will 
designate  desired  speed  as  follov/s: 

1  for  9  niph 

2  for  12  m^ph 

3  for  15  Q'^P'^ 
if  for  13  mph 

Card  form^at  for  these  elements  is  12, 
IX,  12. 

6.  X,  Y  coordinates  for  each  initial  red 

location.   One  card  is  needed  for  each 
red  unit  with  format  F6.1,  IX,  F6.1. 

?•*  If  the  user  has  selected  to  enter  his 

own  routes  in  data  element  5>  the 
following  is  required  for  each  route. 

a.  The  number  of  nodes  in  the  route 
(from  01  to  10).  This  card's  format 
is  12. 
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b,   X,  Y  coordinates  of  each  node 
along  the  route  in  order  from  west 
to  east,   Cne  card  for  each  node  in 
format  F6.1,  IX,  F6.1. 

3.  X,  Y  coordinates  for  blue  location,  force 

level  for  each  unit,  primary  search  direction 
for  that  unit  and  desired  search  width. 
Force  level  must  be  between  1,0  and  3«0j 
search  direction  between  135  and  225  , 
and  search  width  from  30°  to  120  .   One 
card  for  each  blue  unit  with  format  F6.1, 
IX,  F6.1,  IX,  F3.1,  IX,  13,  IX,  13. 

9*  Specify  if  blue  is  to  move  to  alternate 

defensive  positions,  breakpoint  distance, 
and  number  of  10  second  time  intervals 
allowed  for  the  move,   A  1  for  no  move 
or  0  for  option  to  move.   Format  for  this 
card  is  II,  IX,  F6.1,  IX,  12. 

10,"*  Alternate  blue  position  X  and  Y  coordinates 

if  move  is  desired  (data  element  9).  One 
card  for  each  location  with  format  F6.1, 
IX,  F6.1 

11.  Red  weapon  system  hit  and  kill  probab- 

ilities.  The  probabilities  are  entered 
for  6  range  bands,  with  one  card  for  each 

range  band.  Four  probabilities  for  each 

st 
range  are  probability  of  a  1   round  hit. 
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probability  of  a  hit  given  a  previous  round 
hit,  probability  of  a  hit  given  a  previous 
round  miss,  and  probability  of  a  kill 
given  a  hit.  The  6  range  bands  are  from 
0  to  500,  500  to  1000,  1000  to  I5CO, 
1500  to  2000,  2000  to  2500  and  2500  to  3OOO. 
Each  of  these  six  cards  has  the  format 
/+(FZf.2,  IX). 
12.  Same  as  item  11,  except  the  probabilities 

for  the  blue  v/eapon  systems  are  entered. 
Six  cards  with  format  4(^4.2,  IX), 
In  the  formats  listed  above,  a  IX  is  a  space,  an  Ix  is  an 
integer  of  x  digits,  and  a  Fx.y  is  a  real  number  of  length 
X  with  y  digits  to  the  right  of  the  decimal.   For  example, 
an  12  could  be  a  25,  an  F6.1  could  be  3Zf87.4.   A  format 
containing  these  specifications  could  be  12,  IX,  F6.1. 
The  data  card  v;ould  then  look  like  25  3if87.4,  v/ith  the  25 
beginning  in  card  column  one. 
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card  1         2         3         k  iflMent 

column  1234567890123456789012345678901234567890     number 


1  28943 

1 

03  03 

2 

0000.0  2500.0  0500.0  4000.0 

3 

3.0   3-0  3.0 

4 

1  2 

5 

2500.0  7300.0 
2400.0  5900.0 
2450.0  5300.0 


02 

4900.0  7200.0 

7100.0  6700.0 

02 

33^0.0   6000.0  7 

7000.0  6300.0 

01 

6900.0  5500.0 

5700.0  7800.0  3,0   190  120 

5800.0  6500.0  3.0  180  120  8 

5900.0  5900.0  3.0  180  120 


0  0500.0  4  9 


7000.0  7200.0 

7100.0  6500.0  10 

7100.0  5800.0 

0.85  0.S3   0.75  0.70 

0.80  0.80  0.75  0.70 

0.75  0.75  0.70  0.65  11 

0.60  0.65  0.60  0.55 

0.45  0.50  0.50  0.35 

0.20  0.20  0.20  0.20 


Figure  3.   Sample  Data  Deck 


43 


data 
T  -,  ,  element 

card  1  t)  h  nunhpr 

column       123456789012345678901234567890123^567890  nuiaoer 

0.60  0.70  0.65  0.85 

0.85  0.90  0.85  0.90 

0.80  0.85  0.85  0.80  ^^ 

0.75  0.80  0.75  0.70 

0.60  0.70  0.65  0.65 

0.40  0.45   0.40  0.50 


Figure  3»   Sample  Data  Deck  (cont) 
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V.   REQUIRED  JCL 

The  model  and  line-of-sight  data  are  residing  on 
permanent  disk  in  the  '.V.  R,  Church  Computer  Center,   In 
order  to  exercise  the  model,  the  required  job  control 
language  (JCL)  is  illustrated  in  Figure  4.   Due  to  the 
CPU  time  required  for  execution,  a  time  parameter  of 
TIME  =  2  must  be  used  on  the  job  card. 

//JOE  CARD  ,  ti::e=2 

//  EXEC   FORTCLG 

//FORT.SYSIN   DD  UNIT=231if ,  V0L=SER=PAT002,DISP=SHR  , 

//   DSN=S1360.SM0LER.PGM,DCB=(RECFM=FB,IRECL=80,BLKSIZE=3200) 

/* 

//G0.FT08F001  DD  UNIT=2314,V0L=SER=PAT002,DISP=SHR, 

//  DSN=S1360.LOSDATA,DCB=(RECFM=FB,LRECI^80,BLKSIZE=3200) 

/* 

//GO.SYSIN  DD  ^ 

USER  DATA  DECK 

/* 


Figure  k»     Required  JCL 
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VI.   EXPECTED  OUTPUT 

Once  the  model  has  been  executed,  the  user  can  expect 
the  following  output: 

1  -  Program  listing  of  the  model. 

2  -  A  summary  of  the  initial  battle  conditions, 

including  starting  locations  and  options  selected, 

3  -  A  sum.mary  of  battle  conditions  after  each  10 

second  time  step,  including  unit  locations,  force 
levels,  unit  status,  percentage  lost  and  targets 
on  each  unit's  target  list.  The  status  will 
show  one  of  the  following: 

0  -  Unit  alive,  not  firing 

1  -  Unit  alive  and  firing 

2  -  Unit  killed 

3  -  Unit  moving 

4  -  A  statement  of  the  reason  for  battle  termination. 
A  sample  of  the  initial  battle  sum.mary  is  at  Figure  5  and 
the  time  step  summary  is  at  Figure  6. 
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INITIAL  BATTLE   INFORMATION 


UNIT 

X 

Y 

FORCE  LEVEL 

1 

1900.0 

7800.0 

6.0 

2 

1700.0 

6500.0 

6.0 

5 

1700.0 

ii900.0 

6.0 

k 

5800.0 

7700.0 

3.0 

5 

7400.0 

6300.0 

3.0 

6 

6200.0 

5150.0 

3.0 

ATTRITION   IS   STOCHASTIC 

ROUTES   DETERMINED   BY   USER 

RED  VEHICLE   SPEZB  IS   I5.O 

BREAKPOINT   DISTAiNCE   IS   I5OO.O 

BLUE   WILL  MOVE  TO   ALTERNATE   POSITIONS 
ALTERNATE   POSITIONS   ARE: 
UNIT  X  Y 

1+  9300.0        7300.0 

5  9800.0        6200.0 

6  9200.0        5800.0 


RED 

KILL 

PROBABILITIES 

RANGE 

P 

PHH 

PHM 

PKH 

500 

0.95 

0.95 

0.95 

1.00 

1000 

0.80 

0.83 

0.81 

0.95 

1500 

0.65 

0.70 

0.70 

0.90 

2000 

Q.kO 

0.50 

0.55 

0.75 

2500 

0.30 

0.40 

0.40 

0.75 

3000 

0.20 

0.30 

0.25 

0.60 

BLU; 

S  KILL 

1  PROBABILITIES 

RANGE 

P 

PHH 

PHM 

PKH 

500 

0.80 

0.80 

0.80 

1.00 

1000 

0.80 

0.80 

0.80 

0.95 

1300 

0.70 

0.70 

0.70 

0.90 

2C00 

0.65 

0.65 

0.65 

0.85 

2500 

0.45 

0.55 

0.45 

0.70 

3000 

O.i+5 

0.55 

0.45 

0.60 

Figure 

5.   I 

nitial 

Battl 

e  Sum: 
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TIME  =  kSO   SECONDS 


IT 

X 

Y 

FORCE    lEVEL 

STATUS 

LOST-PCT 

TARGETS 

1 

^933.2 

7193.2 

5.0 

0 

0.160 

2 

if832.3 

5799*9 

4.7 

0 

0.214 

3 

4902.6 

4579.8 

6.0 

0 

0.0 

k 

5800.0 

7700.0 

0.0 

2 

1.000 

5 

7400.0 

6300.0 

3.0 

1 

0,0 

2 

6 

6200.0 

5150.0 

3.0 

3 

0.0 

Figure  6.   Time  Step  Summary 
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APPENDIX  3 
Program  Listing 


1  CflMMON  /GflPl/  1PRDIR(6]  ,1SECHD(6)  ,MVTDIfl(6)  ,X(6I  ,T(S)  .SPD(6J 

2  COMMON  /GRP2/  Tfl  (2)  .  Tl  [2)  .  TH  (2)  ,  TM  (2)  ,  TFl  (2)  .  TF2  (3)  ,  TF3  (2)  , 
9           IP  (2,  6)  ,PHH(2,6)  .PHM(2.6)  ,PKH(2,6)  ,TF  (2) 

V  COMMON  /GRP3/  NBU.NRU.FL  (6) , FO  (6) . NOI  (3) , XIC  (3,200) . TIC (3, 200) , 

5  IIOIR (3.200) .RVSP.ISPO 

6  l.lUSTflT (6)  .II (6) .LOST (6.6) . ViSFRfl. VISFRB. SIZETK, 

7  ISIZETW.NT (6) .NF (6) .SRF.DlSMflX. 

8  IMLOSC  16.6)  ,VISFR  (6,6) , RMlNTK.RMXTK.RMlNTM, RMXTW. OP.TOHFR, TNKFR, 

9  IPTTO.S)  .RF.P0fl(6.6)  .ftP0fl(6.6)  .L0fl(6,6)  ,Nfl(6)  .0FL(6)  .POL  (6) 

10  COMMON  /GRPii/  TPOL  (6)  ,aL0O(6.6)  .0(6,6) 

11  COMMON  /GRP5/  LOT  (6. 6) .ROT (6. 6) 

12  COMMON  /HILLS/  XC (100)  . TC (100)  , PERK (100)  . SX (100)  .ST (100)  . RHQ llQO) 

13  COMMON  /HILLS/  SCALE  (100) , TWORHO  (100) , TWOSCL (100) .BASE 
m  COMMON  /HILLS/  NHILLS 

15  COMMON  /COVER/  CXC  (150) , CTC  (150) , CPEflK  (150) . CPXX (ISO) .CPTT  (150) 

IS  COMMON  /COVER/  CPXT  (150) .NCVELS 

17  COMMON  /CflUNTR/  KH. KHH, KV, KN.KGRS, KELL, KINT 

18  COMMON  /GRID/  LST  (10, 10) , NHL  (10. 10) .LISTH  (450) . KHREP (100) .KTREP 

19  COMMON  /GRIQ/  LSTC (10. 10) . NC (10. 10) . LISTC (400)  .KCREP (ISO) 

20  COMMON  /GRP6/  ALPHA  (6) 

21  CflMMON  /GflP7/  XA(6)  ,TA(6)  .IHOVE  (6) 

22  C 

23  C      INITIALIZATION. 

24  C 

25  CALL  OVFLOM 

26  BL-0.0 

27  RL-0.0 

28  HP-0 

29  PflI-22. 0/7.0 

30  ZL». 00001 
91  CALL  SETUP 
32  C 

39  C   READ  TERRAIN  DATA  FOR  LINE  OF  SIGHT 

3«  C   CHECK  FOR  STOCHASTIC  OR  DETERMINISTIC  ATTRITION 

35  C 

36  READ15.13Q)    ITRIT.IS 

37  130  FORMAT (II. IX. IS) 

38  00   132   1-1.6 

39  CALL   RANDOM (IS, TRAN.l) 

«0  flLPHA(l)"{-2.»«TRAN«»<2)  +  (2.  kTRAN*.  3) 

«1  132  CONTINUE 

42  C 

43  C   READ  IN  NUMBER  OF  BLUE  AND  RED  UNITS 

44  C 

45  READ (5.200)  NBU.NRU 

46  200  FORMAT  (12. IX. 12) 

47  C 

48  C   INITIALIZE  WEAPON  SIZES 

49  C 

50  SirETK-2.5 
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51  SIZETH-2.5 

52  C 

53  C   RERO  IN  EFFECTIVE  WEAPON  RANGES 

sv  c 

55  READ(5,102)  RMINTK, RMXTK, RMINTW.RMXTH 

56  102  FORMAT  (F6. 1, IX.F6. I.IX.FS. I.IX.FS.I.IX) 

57  C 

58  C   INITIALIZE  PM, RF, TOWFR. TNKFR  AND  NOD 

59  C 

60  PM=.352 

61  RF-.5 

62  T0WFR=.03 

63  TNKFR».l 
BV  NOD-2 

65  DO  101  I^l.NRU 

66  N0I(n=l25 

67  101  CONTINUE 

68  K-NRU+1 

69  L-NRU+NBU 

70  DO  111  I=1.L  • 

71  ll(l)>0 

72  111  CONTINUE 

73  C 

7H  C   READ  IN  FORCE  LEVELS  OF  EACH  RED  UNIT 

75  C 

76  READ  (5.  103)  (FL  (I)  ,  1  =  1  .NRU) 

77  103  FORMAT  (3  IF3.1, IX)) 

78  C 

79  C   CHECK  FOR  TYPE  OF  ROUTE  DETERHINI TION 

80  C 

81  READ  (5, 106)     IRTE.ISPD 

82  106   FQRMATdl.lX.Il) 

83  IFdSPD.EQ.l)    AVSP=9.0 

84  IF(ISPD.EQ.l)    DST=^0.232 

85  IF (ISP0.EQ.2)    RVSP=12.0 

86  IF  IISPD.EQ.2)    0ST  =  53.6q3 

87  1F(ISP0.EQ.3)    AVSP=»15.0 

88  IFdSPD.EQ.S)    DST=67.053 

89  IF(1SP0.EQ.4)    AVSP=ia.O 

90  IF(ISPD.EQ.ii)    DST=80.il63 

91  C 

92  C   READ  IN  INITIAL  RED  LOCATIONS 

93  C 

9U  DO  6  1  =  1,  NRU 

95  READ(5,107)  XIC  (1 ,  1)  ,  TIC  (1 ,  1) 

96  107  FORMAT  (F6.1, IX, F6.1) 

97  6  CONTINUE 

98  IF (IRTE.EQ.l)  GO  TO  108 

99  DO  2  1  =  1.  NRU 

100  DO  2  J=2,125 


50 


101 

Ticd,  jj-TiCd.j-n-^DSTxCj-n 

102 

xiC(i,j)=xiC(i.j-n-^DSTM(j-n 

103 

IDIR(I,J)-0 

IQil 

2 

CONTINUE 

105 

GO  TO  109 

106 

106 

CALL  ROUTE 

107 

109 

SUMRO»0.0 

108 

00  3  lal.NRU 

109 

FO(l)»FL(n 

110 

SUMR0=SUMR0*F0(1) 

111 

x(n=xic(i.i) 

112 

Y(I)=YIClI.l) 

113 

MVTDIR(I)=IDIR(l,n 

114 

SPDtn-AVSP 

115 

IUSTATm=0 

116 

IPRDIR(I)*IDIR(I.l) 

117 

1SECH0(I)»120 

118 

NFdJ-l 

119 

II  (I]»l 

120 

3 

CONTINUE 

121 

C 

122 

C 

REflO  IN  INITIAL  BLUE  LOCATIONS 

123 

C 

124 

SUMBO=0.0 

125 

00  4  I«K,L 

126 

REA0(5.104)  X  (I)  ,Y  (I)  ,FL  (1)  ,  IPRD 

127 

104 

FORMAT  IF6.1, IX. F6.1,1X.F3.1.1X,1 

128 

FOdJ^FLd) 

129 

SUMB0=SUMBQ-^FOd) 

130 

MVTDIRd)-0 

131 

SPDd)»0.0 

132 

lUSTATd)=0 

133 

IMOVE  d)=0 

134 

4 

CONTINUE 

135 

C 

136 

C 

CHECK  FOR  ALTERNATE  BLUE  POSITIONS 

137 

C 

138 

READ  (5.400]  I  ALT. BREAK. ITEM 

139 

400 

F0RHATdl.lX,F6.1,lX.I2) 

mo 

IF[IALT.EQ.l)  GO  TO  401 

mi 

DO  402  I^K.L 

142 

READ(5,107)  XAd).TA(n 

143 

402 

CONTINUE 

144 

401 

DELTOID. 

145 

TA(1)=20. 

146 

Tl  (1)=8. 

147 

TH(1)«8. 

148 

TMll)*10. 

149 

TFl  (n»l. 

150 

TF2(1)»1. 

ISECUDd) 


AND  READ  IN  IF  WANTED 
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151  TF3(1)-1. 

152  TR(2)-20. 

153  Tl(2)-8. 
ISV  TH(2)>8. 

155  TH(2)-15. 

156  TF1(2)-10. 

157  TF2(2)-12. 

158  TF3(2)-15. 

159  C 

160  C   REflO  IN  HIT  RND  KILL  PROBABILITIES 

161  C 

162  OQ  5  1-1.2 

163  00  514  J»1.6 

1611  REflO(5,5l5]  PlI,  J)  .PHHd.  J)  ,PHM(I,  J)  .PKHd.J) 

165  515  FORMAT  (U  (FU. 2. IX)) 

166  Sm  CONTINUE 

167  5  CONTINUE 

166  PTT(l.n»1.0 

169  PTT(1,2)»0.8 

170  PTT(2.2)=0.2 

171  PTT(1,3)=0.8 

172  PTT  (2,3)  =0.15 

173  PTT(3.3)-0.0T 
17U  DO  31  I-l.NRU 

175  DO  31  J-K,L 

176  NLOSC(I.J)-0 

177  NLOSClJ.I)»Q 

178  Q(I.J1*1.0 

179  QlJ.D'l.O 

180  VISFR(1.J)»0.0 

181  VlSFR(J,I)-0.0 

182  31  CONTINUE 

183  IC-1 
18«    C 

185  C   PRINT  INITIAL  BATTLE  INFORMATION 

186  C 

187  WRITE  (6,599) 

188  599  FORMAT  CI  MX.  'INITIAL  BATTLE  INFORMATION') 

189  WRITE  (6,600) 

190  600  FORMAT  (/IX, 'UNIT'. 7X. "X'.SX. 'T'.qx. 'FORCE  LEVEL*) 

191  DO  601  1=1. L 

192  WRITE(6,602)  I .  X  (I)  ,  T  (I)  ,FL  (I) 

193  602  F0RMflT(lX,I3,3X.F7.1.2X,F7.1,7X.F3.1) 
19U  601  CONTINUE 

195  IF(ITRIT.EQ.l)  GO  TO  603 

196  WRITE(6,60ii) 

197  604  FORMAT  (/IX. 'ATTRITION  IS  STOCHASTIC'/) 

198  CO  TO  605 

199  603  WRITE (6.606) 

200  606  FORMAT (/IX. 'ATTRITION  IS  DETERMINISTIC'/) 
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201  605  IF(IRTE.EQ.O)  GO  TO  607 

202  HRITE (6.6081 

203  608  FORMAT  (IX, "ROUTES  DETERMINED  BY  USER'/) 

204  607  HRITE (6.609)  AVSP 

205  609  FORMAT  (IX. 'RED  VEHICLE  SPEED  IS  '.Fq.l/) 

206  MRITE (6,610)  BREAK 

207  610  FORMAT  (IX. 'BREAKPOINT  DISTANCE  IS  •,F6.1/) 

208  IF(IALT.EQ.O)  GO  TO  615 

209  HRITE  (6,620) 

210  620  FORMAT  (IX, 'BLUE  WILL  NOT  MOVE  TO  ALTERNATE  POSITIONS  V) 

211  GO  TO  625 

212  615  HRITE  (6.630) 

213  630  FORMAT  (IX. 'BLUE  WILL  MOVE  TO  ALTERNATE  POSITIONS '/IX. 
2111  I'ALTERNATE  POSITIONS  ARE: '/IX,  'UNIT '.  5X.  'X  '  ,8X. 'Y') 

215  DO  635  I>K.L 

216  HRITE  (6,6110)  I.XAdl.TAdl 

217  610  FORMAT  (IX, 13, 3X.F7.1.2X.F7.n 

218  635  CONTINUE 

219  625  IRAN=500 

220  HRITE  (6.615) 

221  6U5  FORMAT  (/yx,  'RED  KILL  PROBABILITIES '/IX, 'RANGE '. IX. 'P', 

222  IIX, 'PHH',3X. 'PHM',3X, 'PKH') 

223  DO  650  1-1.6 

2211  HRITE(6.655)  IRAN.P  (1 , 1)  .PHH  (1 ,  I)  ,PHM  (1 ,  I)  ,PKH  (1 ,  11 

225  655  FORMAT  (2X.Ii4,1  (2X.Fq. 2)1 

226  IRAN-IRAN-'500 

227  650  CONTINUE 

228  IRAN>500 

229  HRITE  (6.660) 

230  660  FORMAT  (/4X. 'BLUE  KILL  PROBABILITIES '/IX, 'RANGE '.IX, 'P', 

231  lUX.  'PHH'.SX,  "PMM'.SX,  'PKH') 

232  00  665  I>1.6 

233  WRITE  (6.6551     IRAN.P  (2, II ,PHH  (2, II , PHM  (2. II .PKH  (2, 1) 
2311  IRAN-IRAN-'SOO 

235  665  CONTINUE 

236  HRITE  (6.670) 

237  670   FORMAT  CI '.lOX, 'BATTLE    BEGINS'//) 

238  C 

239  C  UPDATE  LOCATION  OF  RED  UNITS. 
2U0     C 

2m  0ISMAX»5000.0 

2U2  67  DO  9  I-1,NRU 

2U3  IF  (lUSTAT(I)  .EQ.2)  GOTO  9 

244  IF  (lUSTAT  (I) .EQ.O)  GOTO  76 

245  NF  (I)="NF  (I)-»l 

246  IF  (NF(I)  .LT.NOO)  GOTO  9 

247  NF(I1-1 

248  76  00  11  J  -  1.  NRU 

249  IF  (J  .EQ.  II  GO  TO  11 

250  IF  (IUSTAT(J1  .EQ.  21  GO  TO  11 
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251  DIST  -  Xd)  -  X(J) 

252  IF  (DIST  .GT.  DlSMfiX)  GO  TO  9 

253  11  CONTINUE 

254  II  (I)  -  II  (I)  *    1 

255  K7»II(I) 

256  X(n-XIC(I.K7) 

257  T(I)»TIC(I.K71 

258  HVTDIR(I)»IDIFI(I.K7) 

259  IPR0IR(I)«I0IR(I.K7) 

260  9  CONTINUE 

261  C 

262  C  LINE— OF-SIGHT  CHECK  BETWEEN  UNITS  fiNO  TARGETS  SELECTION 

263  C 

26>t  DO  17  J-K.L 

265  NT (J) '0 

266  17  CONTINUE 

267  DO  12  I=1.NRU 

268  NT(I)>0 

269  IF  (lUSTAT  (I) .EQ.2)  GOTO  12 

270  DO  16  J>K.L 

271  IFdUSTRT  (J)  .EQ.2.aR.IUSTnT  (J)  .EQ.3)  GO  TO  16 

272  XXUXdl 

273  TTl-Td) 

27M       133  CALL  ELEV  (XXI , TYl .TMflCI) 

275  XX2-X  (J) 

276  TT2»T  (J) 

277  CALL  ELEV(XX2.TT2.TMACJ) 

278  LATQB=1 

279  LBTOA'l 

260  CALL  L0S(XX1.TT1,TMACI,Q.Q,SIZETK.XX2,TY2,TMACJ.0.0,SIZETW. 

261  ILATOB.LBTOA.VISFRfl.VlSFRB) 
282  VISFR(I.J):iVlSFRA 

263  VISFR(J,n=»VISFRB 

2BH  IF(VISFRA.GT.ZL)  GOTO  18 

285  L0STd.J)«O 

266  LOST  (J, I) »0 

287  NLOSC  d,  J)=NLOSC  d.  J)  ♦I 

288  NLOSC  (J.  I)  =NLOSC  d,  J) 

289  GOTO  16 

290  18  LOST  (I, J) =1 

291  LOST(J,n=>l 
282  NL0SC(1,J)>0 
293  NLOSC(J.n>0 

2911  RRNGE=SQRT  (  (X  (I)  -X  (J)  )  »<«2*  (Y  (I)  -Y  (J)  )  >*x2) 

295  IF  (RANGE. LT.RMINTK. OR. RANGE. GT.RMXTK) GOTO  20 

296  IF(Qd,  J)  .EQ.1.0)  GOTO  20 

297  lUSTflT(n=l 

298  NTd)=NT(n-»l 

299  H^NTd) 

300  LOT  (I. H) -J 
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301  ROT f I, M) -RANGE 

302  IF  (H.EQ.l)  GOTO  20 

303  CALL  SORT  (I.M) 

3011       20  IF  (RANGE. LT.RMINTW. OR. RANGE. GT.RMXTW)  GOTO  16 

305  IF  (Q(J,I)  .EQ.l.O)  GOTO  16 

306  IUSTAT(J)-1 

307  NT  (J)»NT  (J)  *1 

308  H»NT  (J) 

309  LaT(J.H)>I 

310  ROT  {J,M)=RANGE 

311  IF  (H.EQ.l)  GOTO  16 

312  CALL  SORT  (J.H) 

313  16  CONTINUE 
3111        12  CONTINUE 

315  00  25  Ul.NRU 

316  IF  (lUSTATd)  .EQ.2)  GOTO  25 

317  IF  (NT  (I)  .NE.OJ  GOTO  25 

318  IUSTAT(I)-0 

319  NF(n»0 

320  25  CONTINUE 

321  DO  79  J'K.L 

322  IF  (lUSTRT  (J)  .EQ.2.0R.IUSTAT(J)  .EQ.3)  GO  TO  79 

323  IF  (NT  (J)  .EQ.O)  IUSTAT(J)«0 
32*1                  79  CONTINUE 

325  C 

326  C  UPDATE  OF  THE  ACCUMULATED  DETECTION  PROBABILITIES. 

327  C 

328  lAA-1 

329  IBBxNRU 

330  ICC-K 

331  IDD-L 

332  FR-TOWFR 

333  OP=»PM 

334  37  00  m  NIAA.IBB 

335  IF  (lUSTATd)  .EQ.2.0R.lUSTAT(n  .EQ.3)  GO  TO  14 

336  DO  19  J=1CC.100 

337  PROP=0.0 

338  IF(IUSTAT(J)  .EQ.2.0R.IUSTAT  (J)  .EQ.3)  GO  TO  19 

339  OLOQd.  J)  =Q(I,  J) 

SUO  IF  (LOST  (I.J)  .EQ.O)  GOTO  15 

3yi  1F(NT(1)  .GT.OJ  GOTO  22 

342  PCTVIS=V1SFR(J.1) 

3U3  CALL   LAMDAd.J.PCTVIS.DETRAT.PSUBK) 

3«1«1  QV=EXP  (-  (PL  (I)  nDETRAThOPmOELThFL  (J)  )  ) 

3«15  IF  (NT  (J)  .GT.O)    GOTO   23 

3>46  0(1.  J)=Q(I,J)mQv 

3«7  GOTO    19 

348  23   QP«(1.0-PSUBK)«<»«(FRH0ELT«FLd)«FL  (J)) 

349  Qd.  J)=Q(I.  J)m(QV*QP-QVmQP) 

350  GOTO    19 
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351 

22 

N5=NT(n 

352 

DO  24  11=1.N5 

353 

K1=LQT  (].  ID 

354 

flNGl=flTflN2  (Y  (Kl)  -T  (1)  .X  (Kl)  -X  in  ) 

355 

flNG2=flTflN2(T(J)-Y(n  ,X(J)-X(n) 

356 

IF  ((flNGl«flNG2)  .GE.O.Q)  GOTO  77 

357 

IF(flNG2.LT.0.0)  GOTO  32 

358 

flNG=2HPftI*flNGl-flNG2 

359 

GOTO  35 

360 

32 

flNG=2xPflI*flNG2-flNGl 

361 

35 

IF  (fiNG.GT.Pfll)  flNG=2»<PflI-flNG 

362 

GOTO  33 

363 

77 

RNG=RBS(RNG2-RNGn 

3614 

33 

flfl=15.0HPflI/l80.0 

365 

IF  (flNG.GT.flfl)  GOTO  24 

366 

PR0P=PR0P-'PTT[I1.N5) 

367 

24 

CONTINUE 

368 

IF  (PROP. EQ. 0.0)  GOTO  34 

369 

IF  (NT  (J)  .GT.OJGOTO  36 

370 

COLL  LflMDfl(I.J.PCTVlS.DETRflT.PSUBK) 

371 

DETRflT=DETRflTMRF 

372 

QV=EXP  (-  (FL  (I)  «PROP»«DETRflTKOELT«FL  (J)  )  ) 

373 

Qd.Jl  =0(1,  J)  mQV 

374 

GOTO  19 

375 

36 

0(1, J)  =0.0 

376 

GOTO  19 

377 

34 

IF  (Iflfl.EQ.l)  GOTO  19 

378 

0(1.  J)  =1.0 

379 

GOTO  19 

380 

15 

IF(NLQSC(I.J).LE.3)  GOTO  19 

381 

Q(I.J)=1.0 

382 

19 

CONTINUE 

383 

14 

CONTINUE 

384 

IF(Iflfl.EQ.K)  GOTO  38 

385 

FR=TNKFR 

386 

lflfl  =  K 

387 

IBB=L 

388 

ICC=1 

389 

IDD=NRU 

390 

0P=1.0 

391 

GOTO  37 

392 

C 

393 

C 

FIRE  flLLOCflTlON. 

394 

C 

■ 

395 

38 

DO  28  1  =  1. L 

396 

28 

NR  (I)=0 

397 

DO  26  1  =  1. L 

398 

IF  (lUSTRT  (I)  .EQ.2.0R.IUSTRT  (11  .EQ.3)  GO 

399 

IF  (NT  (I)  .EQ.Q)  GOTO  26 

400 

DO  27  J=l,3 

GO   TO   26 


1101  flPQfl(I,J)=0.0 

«02  27    CONTINUE 

«03  IF  (NT  (I)  .EQ.n    GOTO   78 

UQlt  IF  (NT  (I)  .EQ.2)    GOTO   29 

UOS  NOT-3 

406  HMl»L13T(I.n 

407  MM2»LaT(I.2) 
UOe  MH3»L0T(I.31 

■109  PROB=<  (1  .0-Q  (I.HHl)  )mQ  (].HH2)  mQ(].MM3) 

mo  flpofld,  n  =flPOfl  d,  d+ptt  (i,  d  kprob 

411  PRQB=Q(I,MM1)  »<  (1 .0-0  (1 .  MM2)  )  »<Q  ( I ,  MM3) 

412  flP0fl(I,2)=flP0fl(I.2)*PTT  (l.DxPROB 

413  PR08=Q  (I.MMl)  «Q  (I.MM2)  »«  (1.0-0  (1.MM3)  ) 

414  flP0fl(I,3)=flPQfl(1.3)*PTT  (1.1)«PR0B 

415  PRQ8=  (1.0-Q(1,MM1J  )  x  (1.0-0  (I.MH2)  )  «0  (I.MM3) 

416  flPOfl(I,n=flPOfl(I.  1)*PTT  (1.2)hPR0B 

417  flP0fl(1.2)  =flP0fl(I,2)-^PTT  (2.2)  »«PROB 

418  PROB^  ( 1.0-0(1,  HMD  )  xQ  (I.MM2)  x  (1  .0-0  (I.MM3)  ) 

419  flPOfl(l.l)=flPOfl(l.l)-^PTT  (1,2)  kPROB 

420  flPOfl(I.3)=flP0fl(l,3)+PTT  (2,2)«PR0B 

421  PR0B  =  Q(1.MM1)  x  (1 .0-0  (1 .  MM2)  )  »<  (1.0-0  (I,MM3)  ) 

422  flP0fl(I,2)=ftP0fl(I,2)*PTT  (1.2)hPR0B 

423  flP0fl(I,3)«flP0fi(I,3)-^PTT  (2.2)«PR0B 

424  PRQB=  (1.0-0(I,MH1)  )  X  (1.0-0  (1,HH2)  )  «  (1 .  0-0  (I ,  MM3)  ) 

425  flPOfl(l.n-flPOfl(I,l)*PTT  (1.3)xPRQB 

426  flPOfl  (1,2)  =flPOfl  (1.2) +PTT  (2.3)xPR0B 

427  flPOfl  (1,3)  =flPOfl  (1.3)  ♦PTT  (3,3)  xPROB 

428  30  00  114  J=1,N0T 

429  KK  =  L0T(1,J) 

430  Nfl  (KK)  =Nfl(KK)*l 

431  ]N=NR(KK) 

432  L0R(KK,1N)=I 

433  POA(KK.IN)=flPOfl(I,  J) 

434  41  CONTINUE 

435  GOTO  26 

436  29  N0T«2 

437  MM1=LQT(I,1) 

438  MM2=L0T(I,2) 

439  PROB=  (1.0-0(1. MMl)  )  xQ(l.MM2) 

440  flPOfl(I,l)=flPOfl(I.l)*PTT  (I.DkPROB 

441  PR0B=Q(1,MM1)  •«  (1.0-0(1, MM2)  ) 

442  flP0fl(I,2)=fiP0ft  (I,2)*PTT  (1.  DxPROB 

443  PROB=  (1.0-0  (I,MM1)  )  X  (1.0-0  (I,MM2)  ) 

444  flPOfld.  1)  =flPOfl(I.  1)  +PTT  (1.2)  xPROB 

445  flP0fl(I,2)=ftP0fl(l,2)+PTT  12.2)  xPROB 

446  GOTO  30 

447  78  N0T=1 

448  MM1=L0T(I,1) 

449  PRa8=l. 0-0(1. MMl) 

450  flPOfl(l.l)  =flPOfl(l.l)*PTT  (1,1)  xPROB 
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1151 

GOTO  30 

M52 

26 

CONTINUE 

^53 

C 

H5H 

C 

ATTRITION. 

U55 

C 

U56 

SUMR=0.0 

157 

SUMB=0.0 

458 

DO  10  1=1, L 

>i59 

IFClUSTflTd)  .EQ.2.8R.lUSTflT(n  .EQ.3)  GO  TO  10 

ueo 

M6=Nfl(l) 

■161 

SUM=0.0 

162 

IF  (M6.EQ.0)  GOTO  47 

163 

DO  11  J=1,M6 

161 

H7=L0fl(I.  J) 

165 

IF(M7.LT.K)  GOTO  12 

166 

ITTPE=2 

167 

GOTO  13 

168 

12 

ITYPE=1 

169 

13 

RflNGE=>SQRT  (  (X  (I)  -X  IM7)  )  «>«2-f  (T  CI)  -T  (M7)  )  x>*2) 

170 

IF  (ITRIT.EQ.l)  GO  TO  131 

171 

CALL  STQCHnTTPE, RANGE, flJl) 

172 

GO  TO  5000 

173 

131 

CALL  ETKdTTPE, RANGE. T) 

171 

AJI=1.0/T 

175 

5000 

SUH  =  SUM+AJ1><FL  (M7)  xPOA  (I,  J)  xQELT 

176 

11 

CONTINUE 

177 

17 

OFL(I)=FL(I) 

178 

FL  (n=FL(I)-SUM 

179 

IF  (FL  (1)  .GT.ZL)  GOTO  16 

180 

FL  ll)=0.0 

181 

IUSTAT(1)=2 

182 

16 

IF(I.LT.K)  GOTO  60 

183 

SUMB=SUMB*FL  (I) 

181 

TPQL(I)  =  lFO(I)-FLn))/FO(I) 

185 

GO  TO  10 

186 

60 

SUMR=SUMR*FL (I) 

187 

TP0L(n  =  (F0(I)-FL  (DJ/FOd) 

188 

10 

CONTINUE 

189 

C 

190 

C 

PRINT  AND  CHECK  FOR  BATTLE  DETERMINATION. 

191 

C 

192 

ITIME=1Ch10 

193 

DO  57  1=K.L 

191 

IF  IIUSTAT(I)  .EQ.2)  GO  TO  57 

195 

DO  58  J=1,NRU 

196 

IFdUSTAT  (J)  .EQ.2)  GO  TO  58 

197 

CHECK=X  (I)  -X  (J) 

198 

AVD=SQRT  (  (X  (I)  -X  (J)  )  hh2*  (T  (I)-T  (J)  )  m>«2) 

199 

IF  (AVD.LT.BREAK.OR.CHECK.lt. 50.)  GO  TO  250 

500 

58 

CONTINUE 

58 


501  57  CONTINUE 

502  GO  TO  99 

503  C 

504  C        COMPLETE   BLUE   MOVE 

505  C 

506  250  DO   251    l^K.L 

507  IFdflLT.EQ.l.OR.IMOVE  (I)  .EQ.ITEM)    GO   TO   6000 

508  IF  llUSTflT  (I)  .EQ.O)     lUSTflT(l)=3 

509  IMOVE  (I)  >1H0VE  (I)-»l 

510  IFdMOVE  (1)  .LT.ITEM)    GO   TO   251 

511  X(I)=Xfl(n 

512  Y(n=Tfl(l) 

513  IFdUSTflT  11)  .EQ.3)     IUSTflT(n=0 
51V  251    CONTINUE 

515  99   WRITE  (6.112)     ITIME 

516  112   FORMAT  (////IX,  'TIME    =    MU,  IX. 'SECONDS '//) 

517  HRITE(6.113) 

518  113  FORMflTdX. 'UNIT', 5X. 'X', ex. 'T'.SX. 'FORCE  LEVEL  '  ,2X. 'STATUS '. 

519  12X, 'L0ST-PCT',2X, 'TARGETS') 

520  00  59  1  =  1. L 

521  N6>NT(I) 

522  IF  (N6.NE.0)  GO  TO  48 

523  HRITE(6,264)  I .  X  (1)  .  T  (I )  ,  FL  (1)  ,  lUSTAT  (I )  ,TPOL  (1) 

524  264  FORMAT  (3X.I1.3X.F7.1.2X.F7.1.6X.F3. 1.9X.I1.6X,F5.3) 

525  GO  TO  59 

526  48  WRITE  (6.114)  I .  X  (1)  ,  T  (I)  ,FL  (I)  .  lUSTAT  (I)  .TPOL  (I)  , 

527  1  (LOT  (1,  J)  .J=1.N6) 

528  114  FORMAT (3X, 1 1 . 3X.F7. 1 . 2X. F7. 1 . 6X. F3. 1 . 9X. 1 1 . 6X. F5. 3. 3X. 3 (II , IXJ ) 

529  59  CONTINUE 

530  C 

531  C  CHECK  FOR  BATTLE  DETERMINATION. 

532  C 

533  lOT^O 

534  00  53  lal.NRU 

535  IF  (FLd)  .EQ.0.0)  GOTO  53 

536  I0T>1 

537  53  CONTINUE 

538  IFdOT.EQ.l)  GOTO  54 

539  WRITE (6. 117) 

540  117  FORMAT  (IX,  •«».«  RED  FORCE  IS  ELIMINATED.   END  OF  BATTLE.') 

541  GOTO  66 

542  54  lOT^O 

543  DO  55  I=K.L 

544  IF(FLd)  .EQ.O.O)  GOTO  55 

545  lOT^l 

546  55  CONTINUE 

547  IFdOT.EQ.l)  GOTO  65 

548  WRITE  (6. 118) 

549  118  FORMAT  (IX. 'xM«   BLUE  FORCE  IS  ELIMINATED.  END  OF  BATTLE.') 

550  GOTO  66 
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551  6000  WRITE  (6. 119) 

552  119  FQRMflT  (IX,  •«•<«  DISTANCE  BETWEEN  FORCES  IS  TOO  CLOSE.  END  OF  BATTLE 

553  1') 

5514  GOTO  66 

555  65  1C=1C*1 

556  GOTO  67 

557  66  STOP 

558  END 
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1  SUBROUTINE   LOS (Xfl, Tfl. TMflCfl. TMlCfl, SIZEfl, XB, YB.THflCB.TMlCB. SIZES. 

2  -LflTOB.LBTOR.VlSFRfl.VlSFRB) 

3  C 

«  COMMON  /HILLS/  XC  (100)  .  TC  (100)  .PEAK  (100)  .  SX  (100)  .ST  (100)  .RHOIlOOl 

5  COMMON  /HILLS/  SCRLE  (100) , THQRHO (100) , TWQSCL (100] ,BBS€ 

6  COMMON  /HILLS/  NHILLS 

7  COMMON  /COVER/  CXC  (150)  ,  CYC  (150)  .  CPEflK  (150)  ,  CPXX  (150)  .CPTT  (150) 

8  .  COMMON  /COVER/  CPXT  (150) , NCVELS 

9  COMMON  /COUNTR/KH.  KHH.  KV.  KN.  KGRS,  KELL.  KINT 

10  COMMON  /GRID/  LST  (10. 10)  .  NHL  (10. 10)  .  LI  STH  (ti50)  .KHREP  (100)  .KTBCf 

11  COMMON  /GRID/  LSTC (10. 10) . NC (10. 10) . LISTC (400) .KCREP (150) 

12  DIMENSION  IGX  (100)  .  IGY  (100)  .  lEL  (100)  .CSl  (100)  .CS2  UOOl 
19  MTR  NGRIO/10/.GSIZE/lOOO./ 

m  C  SUBROUTINE  TO  COMPUTE  FRflCTIQN  VISIBLE  FOR  OBSERVER  1M\GE.1   PfllW 

15  VISFRA-1. 

16  VISFRB-1. 

17  XBR*X8-XR 

18  YBR-TB-TR 

19  IF((XBR.EQ.O.) .RND.  (YBR.EQ.O.) )  RETURN 
2p  IF(SIZER>TMICR.LE.0.)  GO  TO  510 

21  .  IFISI2EB*TMICB.LE.0.)  GO  TO  510 

22  IFtTHlCfl.LT.O.)  ViSFRfl-l.O+TMlCfl/SIZEfl 
29  IFJTMICB.LT.O.)  VlSFRB-1 .0*TMICB/SIZ£B 
2«  Zfl-TMflCfl  ♦   TMlCfl  ♦   SIZEfl 

25  2B-TMflCB  ♦   TMICB   ♦   SIZEB 

26  KTREP-KTREP*! 

27  ZBR-ZB-ZR 

28  XBRSQ-XBRxXBR 

29  TBRSQ-TBRmTBR 

90  XYBR-XBRmTBR 

91  TU0X8R-2.MX6R 

92  TU0YBR-2.HYBR 

99  C  COMPUTE  GRID  SQUARES  CROSSED  BY  R  TO  B  LINE 

9«l  NGRSQ-0 

95  IF(XBA)  110.95.100 

96  95    XBR-0. 1 

97  100   ISCX— 1 

98  XINC'GSIZE/XBR 

99  GO  TO  120 
«0  110   ISCX-l 

VI  XINC— GSIZE/XBR 

«2  120        1F(T6R)       1140.125.130 

«9  125        TBR-0.1 

«V  130        ISGT— 1 

45  TlNC-GSIZE/TBfl 

46  GO  TO    150 
«7  lUO        ISGT-1 

118  TINC— GSIZE/YBR 

«9  150        IX-1«IFIX(XB/GSIZE1 

SO  IFdX.GT.NGRIOl     IX^NGRIO 
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51  IY-1*1FIX(TB/GS1ZE) 

52  IFCIY.GT.NGRIDl     lY^NGRID 

53  XNEXT«GSIZE«<  iFLOfiT  (IX)  ♦0.5»«  (lSGX-1.)  ) 
SV  TNEXT-GSIZEx  (FLOAT (1 Y) +0. S- (IS6T-1 .)  ) 

55  XSTEP=  (XB-XNEXT) /XB(^ 

56  YSTEP- (YB-YNEXT)/YBfl 

57  160   NGRSQ=NGRSQ-*1 

58  IGX(NGRSQ) -IX 

59  IGT(NGRSQ)-IY 

SO  IF[(XSTEP.G7.1.) .AND.  (YSTEP. GT.l.))  GO  TO  200 

61  IF(XSTEP-YSTEP)  170,180,190 

62  170    IX^IX  -^ISGX 

83  XSTEP-XSTEP-^XINC 

611  GO  TO   160 

65  180   IX-IX-^ISGX 

66  XSTEP=XSTEP*X1NC 

67  190    lY-IY+lSGY 

68  YSTEP-YSTEP-^YINC 

69  GO  TO  160 

70  200   KGRS-K6RS+NGRSQ 

71  C  GRID  SQUARE  LIST  NOU  COMPLETE  IN  IGX.  IGY  UITH  NCRSQ  ENTRIES 

72  C 

73  C  NOW  FIND  WHICH  COVER  ELLIPSES  TOUCH  THE  ft  TO  B  LIHE, 
7«  C  CHECK  ELEVATIONS  AT  SI  AND  S2  FOR  EACH  SUCH  ELLIPSE 

75  NELS-0 

76  CHTMAX-0. 

77  IF(NCVELS.EO.O)    GOTO   270 

78  DO   260   K<=1.NGRSQ 

79  1X»IGX(K) 

80  lY-IGY(K) 

81  N-NC(IX.IY) 

82  IF(N.EQ.O)    GO   TO   260 

83  LS»L3TC(IX,IY) 
8U  LEND-LS+N-1 

85  DO   250   L-LS.LENO 

86  KELL=KELL*1 

87  1C=-LISTC(L) 

88  IF(KCREP(IC)  .EQ.KTREP)    GO   TO   250 

89  KCREP{IC)=KTREP 

90  RX=XA-CXC(IC) 

91  RY-YA-CYC (IC) 

92  PXX=CPXX(IC) 

93  PYY-CPYY(IC) 
9U  PXY=CPXY(IC) 

95  flfl-PXXMXBASQ+PYYxYBASQ+PXYHXYBA 

96  BB-PXXmTWOXBAhRX+PYTmTWOYBAmRY+PXYk (RXKYBfl+RYMXBfi) 

97  CC-PXXmRXxRX+PYYxRYmRY+PXYkRX^RY-1 .0 

98  ARG=BBm6B-U.0mAAmCC 

99  IF(ARG.LE.O.)    GO   TO   250 

100  SQ-SQRT  (ARC) 
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101  SI— (BB-^SQ)/ (S.Owflfl) 

102  S2-(SQ-BB)/(2.0*«RR) 

103  IFtSl.GE.l.)  GO  TO  250 
lOV  IF(S2.LE.O.)  GO  TO  250 
IQS  IFISI.LE.O.)   GO  TO  510 

106  IF(S2.GE.l.)  GO  TO  510 

107  C  NOH  CHECK  LOS  RT  SI  RNO  S2 

108  KINT-KINT*! 

109  CPK-CPEflK(IC) 

110  XS-XR-»S2«<XBR 

111  TS-TR-»S2»«TBfl 

112  CALL  ELEV(XS.YS.HTS) 

113  HTS-HTS-»CPK 
im  ZS-Zfl-S2NZBR 

115  IF(LRTOB.EQ.O)  GO  TO  210 

116  CALL  KQVER(Zfl,TMRCB,SIZEB.ZB.S2.HTS.ZS,VlSFRB) 

117  IF(VISFRB.LE.0.)  GO  TO  510 

118  210   IFILBTQR.EQ.O)  GO  TO  220 

119  S-1.-S2 

120  CRLL  KOVERCZB.TMRCfl.SlZER.Zfl.S.HTS.ZS.VlSFRfl) 
IZU.  IF  (VISFRfl.LE.O.)  GO  TO  510 

122  220   XS-XR-^SImXBR 

123  YS-TR-»S1kTBR 

1211  CALL   ELEV(XS.YS.HTSJ 

125  HTS-HTS+CPK 

126  ZS-ZR-^SImZBR 

127  IF  (LRTOB.EQ.O)  GO  TO  230 

128  CALL  KQVER(Zfl.TMRCB.SIZEB.ZB.Sl,HTS,ZS.VlSFRB) 

129  1F(V1SFRB.LE.0.)  GO  TO  510 
13Q  230   IFtLBTOR.EQ.O)  GO  TO  240 

131  S>1.0-S1 

132  CALL  KOVERlZB.TMRCfl.SIZEA.ZA.S.HTS.ZS.VlSFRA) 

133  IF  (VISFRfl.LE.O.)  GO  TO  510 
13U  2«0   NELS-NELS^l 

135  IEL(NELS)>IC 

136  CSllNELSl^Sl 

137  CS2(NELS)-S2 

138  IF  (CPK.GT.CHTMflX)  CHTMflX  =  CPK 

139  250   CONTINUE 

140  260   CONTINUE 

lUl  C  ALL  ELLIPSES  CHECKED 

1>12  C 

m3  C  NOW  START  ON  THE  HlLLS 
IVV      270  00  600  K^l.NGRSQ 

145  IX'IGX(K) 

146  IT=IGY(K) 

147  IF(NHL(IX,IY1  .EQ.O)  GO  TO  600 

148  LS»LSTIIX.IY) 

149  LEND=LS*NHL  (IX.IY)  -1 

150  DO  500  L»LS,LEN0 
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151  I-LISTH(L) 

152  IFCKHREPd)  .EQ.KTREP)    G8   TO   500 

153  KHREP(I)-KTREP 

1S«    C  PROCESSING  FOR  HILL  I  STARTS  HERE 

155  KH"KH-»1 

156  C  COMPUTE  H  -TOP  OF  THIS  HILL  ALONG  0-T  LINE 

157  CX-XBA/SX(1J 
15a  CY-TBA/ST(I) 

159  OX- (XA-XC(I)  ) /SX(I) 

160  DT- (YA-YCd)) /STd) 

161  FQ-THOSCL  U)  «  (CXxDX-^CYhQT+RHO  (I)  •«  (CXnOY+CThDW  1 

162  GQ=SCflLE  (I) « (CXxCX+CYmCY+TWQRHO (I)  xCXmCT) 

163  1F(GQ  .EQ.O.)   GO  TO  500 
16M  N— FQ/(2.kGQ) 

165  IF(ABS(H)  .GT.5.)  GO  TO  500 

166  FSQ^FQmFQ 

167  EC-SCALE  (1)  »«  (OXmOX+OYxDY  +  TWORHO  ll  1  «DX»<OT) 
188  PQHER-EQ-FSQ/(4.mGQ) 

163  IF  (POWER   .LT.    -3.)    GO   TO   500 

170  HHH-PEAK  in  "EXP (POWER] 

171  KHH-KHH+1 

172  IF(HHW.LE.BASE)  GO  TO  500 

173  ZW-ZA-»HmZBR 

17«l  IF  (  (M.LT.O.)  .OR.  (H.GT.l.)l  GO  TO  300 

175  IF(HHW  .GE.  ZW)  GO  TO  510 

176  CVHTW-O. 

177  IF(NELS.E(3.0)  GO  TO  300 

178  00  280  H-1,NELS 

179  IF  ((CSl  (H) .GE.W) .OR.  (CS2  (H) .LE.M))  GO  TO  260 
laO  IC-I£L(H) 

181  IF(CVHTW.LT.CPERK{IC))   CVHTW-CPEAK  (IC] 

182  280   CONTINUE 

163  IF((HHH*CVHTW) .GE.ZWJ   GO  TO  510 

lan     300   IF(HHW*CHTMAX.LT.RMIN1  lZA-SIZEfl,ZB-SlZEB)l  GO  TO  500 

185  C  IF  HE  GET  TO  HERE  THEN  NEED  TO  FIND  LOWEST  SIGHT  LINE  OVER  HILl 

186  C  NEHTON  ITERATION  A  TO  B  GIVING  VISFRB 

187  IF(LATOB.EQ.O)  GO  TO  UOO 

188  KV-KV+1 

189  V-H 

190  HHV-HHM 

191  NCT-0 

192  FV'FQhV 

193  TW0GV-2.mGQmV 

194  330    FCNV-ZR*HHV»«(TW0GV«<V*FV-1.) 

195  KN^KN^l 

196  FflCTORa  (THOGVxTWOGV+2. H (GQ+TWOGV«FQ) ♦FSO) 

197  DFCNV-HHVmV«FACTOR 

198  lF(flBS(DFCNV)  .LT.l.E-10)  GO  TO  350 

199  V-V-FCNV/DFCNV 

200  IF(RBS(V)  .GT.S.IGO  TO  >400 
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201  FV-FQmV 

202  TU0GV-2.»GQmV 

203  POWER  -  EQ*FV*GQ«V««V 

204  IF  (POWER  .LT.  -3.)  GO  TO  UOO 

205  HHV-PEflKlDxEXPCPOWER) 

206  DHHV-HHVw  (FQ+TWQGV) 

207  ELV-Zfl+DHHVHV 

208  IF  IflBS(HHV-ELV)  .LT.l.)  GO  TO  350 

209  NCT-NCT+l 

210  IF  (NCT.LT.IO)  GO  TO  330 

211  350   IF((V.LT.O.)  .OR.  (V.GT.l.l)  G0Ta>400 

212  CVHTV«0. 

213  IF(NELS.EQ.03  GO  TO  390 

214  00  380  H-l.NELS 

215  IF  (  (CSl  IM) .GE.V) .OR.  (CS2  (H) .LE.V) ) GO  TO  380 

216  IC=IEL(H) 

217  IF  ICVHTV.LT.CPEflK  llO)  CVHTV-CPEflK  (IC) 

218  380   CONTINUE 

219  390    HTV»HHV*CVHTV 

220  ZV=Zfl-»VMZBfl 

221  CALL  KOVER(Zfl,TMflCB,SlZEB,ZB,V.HTV,ZV.VISFRB) 

222  IF  (VISFRB.LE.O.)  GO  TO  510 

223  C  NEHTQN  ITERRTION  B  TO  fl  GIVING  VISFRR 

224  400   IF(LBTOP.EQ.O)  GO  TO  500 

225  KV-KV+1 

226  V>W 

227  VM1»V-1. 

228  HHV-HHW 

229  NCT-0 

230  FV-FQmV 

231  TM0GV:>2.mGQmV 

232  430   FCNV-ZB*MHV»«  (  (FQ+TWOGV)  mVMI-1.1 

233  KN»KN-»1 

234  FflCTOR»  (TW0GV«TW0GV*2.  «  (GQ*TWOGV«FQ)  ■►FSQ) 

235  DFCNV-HHVwVMlxFflCTOR 

236  IFlflBSlDFCNV)     .LT.l.E-10)    GO   TO   450 

237  V-V-FCNV/DFCNV 

238  IFlflBS(V)  .GT.5.)G0  TO  500 

239  VMl-V-1. 

240  FV^FQkV 

241  TM0GV>2.mGQmV 

242  POWER  -  EQ*FV+GQ»«V«V 

243  IF  (POWER  .LT.  -3.)  GO  TO  500 

244  HHV»PEfiK(nHEXP(POHER) 

245  DHHV=HHVk  (FQ+TWOGV) 

246  ELV=ZB-»DHHV*iVMl 

247  IF    (flBSlHHV-ELV)     .LT.l.)    GO   TO   450 

248  NCT>NCT-»1 

249  IF  (NCT.LT.IO)  GO  TO  430 

250  450   IF((V.LT.O.) .OR.  (V.GT.l.) )  GO  TO  500 
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251 

CVHTV-0. 

252 

IF  (NELS.EQ.O)  GO  TO  190 

253 

DO  180  M-»1,NELS 

254 

IF  (  (C31  (M)  .GE.  V)  .OR.  (CS2  (M)  .LE.V)  )  GO  TO  ISO 

255 

IC»IEL(M) 

256 

IF  (CVHTV.LT.CPEflK (IC) )  CVHTV  =  CP£flK (IC) 

257 

U80 

CONTINUE 

258 

490 

HTV=MHV*CVHTV 

259 

ZV-Zfl+VxZBfl 

260 

S»-VM1 

261 

CALL  KOVER(ZB.TMflCfl,SIZEfl,Zfl,S,HTV.ZV.VISFRfl) 

262 

IF  (ViSFRfl.LE.O.)  GO  TO  510 

263 

500 

CONTINUE 

2611 

600 

CONTINUE 

265 

RETURN 

266 

510 

VISFRfl=0. 

267 

VISFRB=0. 

268 

RETURN 

269 

END 
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1  SUBRQUTINE  SETUP 

2  C 

a  COMMON  /HILLS/  XC  (100)  .  TC  (100)  ,  PEAK  (100)  ,  SX  (100)  ,  ST  (100)  .RHO  (100) 

U  COMMON  /HILLS/  SCALE  (lOO)  , TWQRHO  (1 00)  , TWOSCL  (100) , BASE 

5  COMMON  /HILLS/  NHlLLS 

6  COMMON  /COVER/  CXC (150) ,CTC (150) . CPEflK (150) ,CPXX  (150) , CPTY  (150) 

7  COMMON  /COVER/  CPXT  (150) , NCVELS 

8  COMMON  /CQUNTR/KH,KHH,KV,KN,KGRS,KELL,KINT 

9  COMMON  /GRID/  LST  (10, 10) , NHL  (10, 10) . LISTH  (450) ,KHREP  (100)  , KTREP 

10  COMMON  /GRID/  LSTC (10, 10) , NC (10, 10) . LISTC  (UOO) . KCREP  (150) 

11  READ  (8. 7)  NHILLS 

12  READ  (8.47)  BASE 

13  47    FORMAT  (FIO. "4) 
IV  7     FORMAT  (16) 

15  17    F0RMAT(5F8.3,F6.4) 

16  00  50  I-l, NHILLS 

17  READ(8,17)  XC  (I)  ,  TC  (I)  ,  PEAK  (I)  .  SX  (I)  .  ST  (I)  ,  RMO  (I) 

18  50    CONTINUE 

19  READ(8,37)LST 

20  READ  (8. 37) NHL 

21  REA0(8,7)NHT0T 

22  REA0(8.37)  (LISTH  (I)  ,  1  =  1  .NHTQT) 

23  37    FORMAT (1615) 

24  REAQ  (8.71  NCVELS 

25  IF  (NCVELS. EQ. 01  GO  TO  65 

26  DO  60  I-l. NCVELS 

27  REA0(8.27)CXC(I)  ,CTC(I)  .CPEAKd)  .CPXXd)  .CPTT(I)  .CPXT(I) 

28  27    FeRMAT(3Fl0.i4.3El3.7) 

29  KCREP  (I) --21471183600 

30  60    CONTINUE 

31  RER0(8.37)LSTC 

32  READ  (8. 37) NC 

33  RERD(8.7)NCT0T 

3«*  READ(8,37)  (LISTC  (I)  .  I-l  .NCTQT) 

35  85  00  100  I-l. NHILLS 

36  SX(I)-SX(I)Ml.625 

37  ST(I)-ST(I)Ki.625 

38  XC(I)»(XC(I)-500.)«100. 

39  TC(1)-(TC(I)-930.)>.100. 
VO  TWORHO(I)  =2.»«RH0(I) 

m  SCALE(I)— 1./  (2.«(1.-RH0(I)»«k2)  ) 

42  TH0SCL(I)-2.«SCALE(n 

•13  KHREP  (I)— 2147483600 

44  C   ALL  VALUES  NOW  IN  METERS  ON  0  --  10,000  GRID 

45  100   CONTINUE 

46  KTREP— 2147483600 

47  KH-0 

48  KHH-0 

49  KV-0 

50  KN-0 
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51 

KGRS-0 

52 

KELL=Q 

53 

KlNT-0 

54 

RETURN 

55 

END 
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1  c 

2  SUBROUTINE  ROUTE 

3  C 

II  COHHON  /GRP3/  NBU.  NRU.FL  (6)  .FO  (6)  .  NOI  (3)  .  XIC  (3.200)  ,  TIC  (3,200} . 

5  1IDIR(3.2Q0) ,flVSP.ISPD 

6  l.lUSTflT  (6) .11  (61 .LOST  (6.6) . VlSFRfl, VISFRB, SIZETK, 

7  ISIZETW.NT (6) ,NF (6) .SRF.DISMflX. 

8  INLOSC (6.6) .VISFR  (6.6) .RMINTK, RMXTK. RMINTH, RMXTH, OP, TOWFR.TNKFB, 

9  IPTT  (3.3)  .RF.POfl  (6,6)  .flPOfl  (6,6)  .LOfl  (6.6)  ,  Nfl  (6)  .OFL  (6)  .POL  (6) 

10  DIMENSION  XL0C(3.20)  .TL0C(3,20)  .N(3) 

11  IFdSPQ.EQ.^)  0ST>80.>J63 

12  1F(ISPD.EQ.3)  DST  =  67.053 

13  IF(ISPD.EQ.2)  DST=53.643 
1«  IF  (ISPO.EQ.l)  DST.40.232 

15  DO  300  I-i.NRU 

16  RERD(S.15)  N(I) 

17  15  FORMAT (12) 

18  NLoNd)^! 

19  DO  200  IN>2,NL 

20  RERO  (5.201)  XLOCS.TLQCS 

21  201  FORMAT (F6.1. IX. F6.1) 

22  XLOC(I.IN)>XLaCS 

23  YLaC(I,IN)-YLOCS 

24  200  CONTINUE 

25  XLOCd.D-XICd.l) 

26  TLOC(I.l)-TlCd.l) 

27  IOIRd.l)>0 

28  NL-N(I) 

29  NUM«2 

30  DO  305  J-l.NL 

31  XL-XLOCd.J+D-XLOCd.J) 

32  TL»YL0C(1,J*1)-YL0C  (I.J) 

33  D1ST»SQRT  (XLm«2-»YL»«»«2) 

34  Y-flBS(YL) 

35  Z»Y/XL 

36  RNGL>)^TRN(Z) 

37  DEG-flNGLH57.2958 

38  IF  (J.EQ.l)  GO  TO  320 

39  XLN=(DST-EXTRfl)HCOS(fiNGL) 

40  DIST=  (DIST+EXTRfl)-DST 

41  YLN=(DST-EXTRfl)«SlN(flNGL) 

42  XICd.NUM)  »XIC(I.NUM-1)  ♦XLN-^XLE 

43  IF(YL.GT.O.)    GO   TO   325 

44  YLN—YLN 

45  325   YICd.NUM)=iYIC(I.NUM-l)-»YLN*YLE 

46  1F(YL.GT.0.)    GO   TO   340 

47  IDIRd.NUM)— IFIX(DEG) 

48  GO   TO  341 

49  3>10    IDIRd.NUM)  »IFIX(0EG) 

50  341    NUM>NUM-»1 
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51 

320 

XLN=DSTKCOS(flNGL) 

52 

TLN=DST«SlN(flNGL] 

53 

IFITL.GT.O.)  GO  TO  310 

5«t 

TLN—TLN 

55 

310 

IF(DIST.LT.DST)  GO  TO  315 

56 

XIC(I.NUM)»XIC(I,NUM-1)  ♦XLN 

57 

TIC(I.NUM)-TlC(I,NUM-l)*rLN 

58 

IFCYL.GT.O.)  60  TO  342 

59 

ID1R(I.NUH)=.-IFIX(DEG) 

60 

GO  TO  343 

61 

342 

IDIR(I,NUM)»IF1X(DEG) 

62 

313 

DIST=013T-DST 

63 

NUM-NUM+1 

64 

60  TO  310 

65 

315 

EXTRfl>OIST 

66 

XLE=EXTRft«COS(flNGL) 

67 

YLE=EXTRflMSlN(flNGL) 

68 

IFITL.GT.O.)  60  TO  305 

69 

YLE=-TLE 

70 

305 

CONTINUE 

71 

300 

CONTINUE 

72 

RETURN 

73 

END 
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1  c 

2  SUBROUTINE    LflMDfl  (I . J, PCTVIS, DETRflT, PKV 

3  C  / 

V      C  SUBROUTINE  TO  COMPUTE  DETECTION  RATE /(DETRflT)  OF  TARGET  J 

5  C  BT  OBSERVER  I  GIVEN  THE  VISIBLE  FRflCtlON  (PCTVIS) . 

6  C 

7  COMMON   /GRPl/    IPRDIR{6)  .ISECH0(6)  ,MVTDIR(6)  .X(6)  ,Y(6)  ,SPD(6) 

8  TCFflCT-1.0 

9  ZEROL=>0. 00001 

10  Pfll«22. 0/7.0 

11  7      D-nSECWD(l)«Pfll/180. 01/2.0 

12  BBB«  (1.0/(2.0m  (SIN(D)-OMCOS(On)) 

13  IF(flBS(BBB)  .LT.ZEROL)    BBB^O.O 
IV  flAR>  (-B6B)xC0S(0) 

15  IF  (flBS(flflfl) .LT.ZEROL)    flflfl-0.0 

16  aTflNG*flTRN2  (  (T  (J)  -T  (I)  )  ,  (X  (J)  -X  (I)  ) ) 

17  PD»IPRDIR(I)»<PflI/l80.0 

18  IF  ((POxOTRNGl  .GE.0.0)  GOTO  1 

19  IF  (PO.LT.O.O)  GOTO  9 

20  flNGLE-2>«PflI*0TflNG-PD 

21  GOTO    IQ 

22  9   flNGLE=2»«PflI*PD-aTflNG 

23  10    IF  (ANGLE. GT.Pfll)    flNGLE=i2MPAI -ANGLE 
2>l  GOTO  2 

25  1    ANGLE-ABS(PO-OTANG) 

26  2    IF (ANGLE. GT.O)    GOTO   3 

27  OUP-PD+0 

28  DLOH-PO-D 

29  ANGLFT-OTANG+dS.O-PAI/iaO.) 

30  IF(ANGLFT.GT.DUP)    RNGLFT=DUP 

31  ANGLRT=»OTRNG- (15.»«PAI/180.) 

32  IF (ANGLRT.LT.DLOW)  ANGLRT=DLOW 

33  PK-BBBk  (SIN  (ANGLFT) -SIN  (ANGLRT) ) ♦AAA- (ANGLFT-ANGLRT) 
Sq  IF  (PK.LT.O.O)  GOTO  3 

35  IF (PK.GT.l.O)  GOTO  5 

36  GOTO  8 

37  3  PK=0.0 

38  DETRAT=0.0 

39  GOTO  6 
t|0  5  PK-l.O 

111  8   RANGE=SQRT  ( (X  (J)  -X  (I) )  >*h2-^  (T  (J)  -T  (I)  )  •<x2) 

V2  RR-O.OOl-RANGE/PCTVIS 

«3  T0ANG-ATAN2(  (T  (n-T(J))  ,  (X  (n-X(J))  J 

UV  flD=MVTDIR(J)MPAI/l80.0 

«5  HORVEL«ABS (SPD (J) xSIN (TOANG-AD)  ) 

«6  DEN0M=l.il53*TCFACTx  (0.  5978*2.  188»«  (RR»<»«2)  -0.  5038«H0RVEL) 

U7  IF(DENOM.LE.ZEROL)    DENOM-ZEROL 

M  DETRAT=0.003*1.088/DENOM 

•19  DETRAT=DETRAT«PK 

50  6   RETURN 
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1  SUBROUTINE  ELEV  (X, T, TMflC) 

2  C 

3  COMMON  /HILLS/  XC  (100)  .  TC  (100)  .  PEAK  (100)  ,  SX  (100)  ,  ST  (100)  .RMQ  (IQQ) 

4  COMMON  /HILLS/  SCALE  (100) , TWORHO  (100) , TWOSCL (100) , BASE 

5  COMMON  /HILLS/  NHILLS 

6  COMMON  /GRID/  LST  (10, 10)  .  NHL  (10,  10)  ,  LISTH  (450)  .  KHREP  (100)  ,KTREP 

7  COMMON  /GRID/  LSTC  (10, 10) , NC  (10. 10) , LISTC  (400) , KCREP (ISO) 

8  DflTfl  NGR1D/10/,GSIZE/1000./ 

9  C  FUNCTION  TO  COMPUTE  TERRAIN  ELEVATION  FOR  GIVEN  X,  Y  COORDINATES. 

10  ZMAX^BASE 

11  1X-1*IFIX(X/GSIZE) 

12  IF(IX.GT.NGRID)   IX=NGRID 

13  IT-1*IFIX(T/GSIZE) 

m  IF  (lY.GT.NGRID)   IT=NGR1D 

15  IF  (NHL(IX,IT)  .EQ.O)  GO  TO  150 

16  LS=LST(IX.IT) 

17  LENO=LS*NHL (1X,IT)-1 

18  00  100  L'LS.LENO 

19  I-LISTH(L) 

20  QX- (X-XCdl  ) /SX(I) 

21  QXSQ<:QXmQX 

22  IF  (QXSQ  .GE.9.)  GO  TO  100 

23  or- (Y-TCII))/SY(1) 
21  QYSQ-QYwQY 

25  IF  (QYSQ  .GE.  9.)  GO  TO  100 

26  QXY-TWORHO(I)»<QX«QY 

27  FACTaR»SCALE(I)»<(QXSQ+QYSQ-»QXY) 

28  IF  (FACTOR. LT. -3.)  GO  TO  100 

29  HT-PEAK (I) «EXP (FACTOR) 

30  IF  (HT.LE.ZMAX)  GO  TO  100 

31  ZMAX  »  HT 

32  100    CONTINUE 

33  150    TMAC=ZMAX 
3U  RETURN 

35  ENO 
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1  SUBROUTINE    STOCH  (I . RflNGE.fl) 

2  C 

3  C    SUBROUTINE  TO  COMPUTE  STOCHASTIC  ATTRITION  COEFFICIENT 
•1  C 

5  COMMON  /GRP6/  flUPHfl  (6) 

6  COMMON  /GRP3/  NBU.NRU, FL  (6) . FO  16) , NOI  (3)  ,XIC (3,200) , TIC (3,2QQJ , 

7  1ID1R(3.200) .flVSP.lSPO 

8  l.lUSTflT  (6) . II  (6) .LOST (6.6) . ViSFRfl. VISFRB.SIZETK, 

9  1S1ZETH,NT(6)  .NF(6)  .SRF.DlSMflX. 

10  INLOSC  (6.6) . VISFR  (6.6) . RHINTK. RMXTK. RMINTW. RMXTM, OP. TOWFR, TNKFR, 

11  1PTT(3,3)  ,RF,P0fl(6.6)  .flP0fl(6.6)  ,L0fl(6.6)  .Nfl(6)  ,0FL(6)  ,P0L(6J 

12  IF(I.E(3.2)  GO  TO  5003 

13  fl-flLPHfl  (I) X  (  (1 .0-RflNGE/RMXTW) m«2) 
IV  GO   TO   50GU 

15  5003   fl-flLPHfl  (I)  »<  ((1.0-RflNGE/RMXTK)  x»«2) 

16  50014   RETURN 

17  END 
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1  SUBROUTINE  ETK  (I , RANGE. T) 

2  C 

3  C     SUBROUTINE  TO  COMPUTE  EXPECTED  TIME  TO  KILL  fl  TARGET. 
II      C 

5  COMMON  /GRP2/  Tfl  (2)  ,  Tl  (2)  .  TH  (2)  .  TM  (2)  .  TFl(2)  ,  TF2  (2J  ,  TF3  (2)  , 

6  IP  (2,6)  ,PHH12,6)  ,PHM(2.6)  .PKH(2,6)  ,  TF  (2) 

7  IF(I.EQ.2)  GOTO  5 

8  TFin-TFMl) 

9  GOTO  6 

10  5  IF  (RANGE. GT.IOGO.OJ  GOTO  7 

11  TF  (D-TFl  (I)-  (TFl  (IJx  (lOOO.O-RANGE)  /1 000. 01 

12  GOTO  6 

13  7  IF  (RANGE. GT. 2000.0)  GOTO  6 

IV  TF  (I) -TF2  (I) -  (  (TF2 (I) -TFl (1) ) k (2000.0-RflNGE)  /lOOO.O) 

15  GOTO  6 

16  8  TFlI)«TF3(n-((TF3(I)-TF2(I))x(3000.0-RflNGE)/lQ00.G) 

17  6  J-(RRNGE*25Q.O)/500.0 

18  IF(J.GT.6)  J=6 

19  T-TA(1)*T1  (D-THd)  ♦((TH(1)*TF  (I)  )  /PKH(I,  J))  +  (  (TM  (I)  ♦TF  (I) )  / 

20  iPHMd,  J)  )  >»{  (l.Q-PHH  (I.J)  )/PKH(I.  J)  ♦PHHd,  J)  -P  (1,  J)J 

21  RETURN 

22  END 
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1  c 

2  SUBROUTINE  SORT  (I, M) 

3  COMMON  /GRPS/  LOT  16, 6) , ROT  (6, 6) 

4  DO  19  J:>1,M 

5  IF(ROT(I.MJ  .GE.ROTd.JJ)  GOTO  19 

6  21  R-ROT  (I.J) 
7 
8 
9 
10 

n 

12 
13 
lU 


NN=LOT(I.J) 

ROT(I.J)=ROT(I,M) 

LOTlI,J)=LaT(I,M) 

RQT(1,H)«R 

LOT  (I,M)-NN 

19 

CONTINUE 

RETURN 

END 
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1  SUBROUTINE  KOVER(ZO,TMflCT.SlZET,ZT,S.HTS.ZS,VlSFRT) 

2  C 

3  IF(S.EQ.O.)  G8  TO  2000 

«  IF  (HTS.GE.ZS)  GO  TO  2050 

5  HEXT  =  ZO-^(HTS-ZO)/S 

6  EViST-flMflXl  lHEXT,TMflCT) 

7  IF  (EVIST.6E.ZT)  GO  TO  2050 

8  IF  (EVIST.LE.ZT-SIZET)  RETURN 

9  VIS=  (ZT-EV1ST)/SIZET 

10  IF(VIS.LT.VISFRT)  VISFRT=VIS 

11  RETURN 

12  2000   IF (HTS.LT.ZQ)  RETURN 

13  2050  VISFRT=0.0 
IV  RETURN 

IS  END 
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APPENDIX  C 
Definition  of  Variables  in  Computer  Program 


ALPHA(I)     =  Initial  attrition-rate  coefficient  for 

stochastic  attrition  module. 

th 
APOA(I,J)    =  The  average  proportion  of  the  j   attacker 

of  unit  i  allocated  to  fire  on  unit  i, 

AVSP         =  Average  speed  of  moving  attacking  units, 

BREAK        =  Breakpoint  distance  between  attackers  and 

defenders. 

DISMAX       =  Maximum  distance  allowed  between  attacking 

units  before  the  leading  unit  is  delayed. 

DIST         =  The  straight-line  distance  between  two 

movement  nodes  inputed  by  the  user. 

DST         =  The  distance  to  be  moved  each  time  step 

by  attacking  units.  ^ 

FL(I)        =  Force  level  of  unit  i. 

FO(I)        =  Initial  force  level  of  unit  i. 

lALT         =  Denotes  whether  the  user  desires  alternate 

defensive  positions  or  not. 

IC  =  Counts  number  of  time  units  a  defender 

has  been  moving. 

IDIR(I,J)    =  Direction  of  j^^  interval  in  i^^  route. 

II(I)        ■  Interval  index  for  unit  i. 

IMOVE        =  Number  of  time  units  a  defender  is  allowed 

for  moving  to  an  alternate  position. 

IPRDIP(I)    =  Primary  direction  of  movement  for  unit  i. 

IRTE        =  Denotes  whether  user  v/ants  to  input  routes 

or  not. 

IS  =  Random  number  seed  used  for  stochastic  attrition, 

ISECWD(I)    -  V/idth  of  search  sector  for  unit  i. 
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ISPD 

ITEM 

ITIME 
ITRIT 

lUSTAT(I) 

IjCA(I,J) 

LOST(I,J) 

LC'T(I,J) 

MVTDIR(I) 

N(I) 

NA(I) 

NBU 
IIF(I) 

NLOSC(I,J) 

NOI(I) 
IJRJ] 
NT(I) 
OFL(I) 

P(I,J) 
PHH(I,J) 


=  Input  variable  to  denote  user's  desired 
speed  for  attackers  movement. 

=  Input  variable  denoting  number  of  time 
steps  allowed  for  defender's  move. 

=  Current  time,  in  seconds,  of  battle. 

=  Input  variable  denoting  whether  attrition 
will  be  stochastic  or  deterministic. 

=  Current  status  of  unit  i. 

th 
=  The  number  of  the  j   attacker  of  unit  i. 

=  Denotes  whether  line-of-sight  exists 
between  unit  i  and  3  or  not. 

=  The  number  of  the  j   target  of  unit  i. 

=  Movement  direction  of  unit  i. 

=  Number  of  nodes  inputed  by  user  for 
route  i. 

=  Number  of  attackers  of  unit  i, 

=  Number  of  blue  units. 

=  Number  of  time  units  unit  i  is  allowed  to 
fire  at  the  same  location. 


line-of-sight  does  not  exist  between  unit 


th 


Num.ber   of   continous   time   stei3s   that 

line-of-sight 

i  and  unit  j. 

Number  of  intervals  in  the  i"'^  route. 

Number  of  Red  Units. 

Number  of  targets  of  unit  i. 

Force  level  of  unit  i  during  previous 
time  step. 

st 
Probability  of  1   round  hit  by  unit  i 

in  range  band  j. 

Probability  of  a  hit  following  a  hit 
by  unit  i  in  range  band  j. 
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PHM(I,J)     =  Probability  of  a  hit  following  a  miss 

b"  unit  i  in  range  band  j, 

PimdjJ)     =  Probability  of  a  kill  given  a  hit  by 

unit  i  in  range  band  j, 

PM  =  The  proportion  of  time  a  moving  unit 

is  searching  for  targets, 

POA(I,J)     =  The  proportion  of  the  j  ^  attacker  of 

unit  i  allocated  to  fire  on  unit  i. 

POL(I)       •  Percent  of  unit  i  lost  during  the 

current  time  step, 

?TT(I,J)     =  Proportion  of  surviving  firepower 

allocated  to  the  i''^^'^  target  if  there 
are  j  targets  available. 

RANGE        =  Current  minimum  distance  between  attackers 

and  defenders. 

RMINTK       =  Minimum  effective  range  for  attacking 

weapon  system., 

RMIITT'.V       =  Minimum  effective  range  for  defending 

weapon  system. 

RMXTK        =  Maximum  effective  range  for  attacking 

v/eapon  system. 

RMXTW        =  Maxim-um  effective  range  for  defending 

weapon  system. 

RCT(I,J)     a  The  range  of  the  j^^  target  of  unit  i, 

SIZETK       =  Size  of  attacking  vehicle. 

SIZETV/       =  Size  of  defending  vehicle. 

TA(K)        =  Time  to  acquire  a  target  for  k  '  weapon 

system  type  (k  =  1,2). 

TFl(K)       =  Tim.e  of  flight  to  1000m  for  k^^  weapon 

system  type  (k  =  1,2). 

TF2(K)       =  Time  of  flight  to  2000m  for  k^^  weapon 

system  type  (k  =  1,2). 

TF3(K)       =  Time  of  flight  to  3000m  for  k^^  weapon 

system  type  (k  =  1,2). 
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TH(K) 

TI(K) 

TM(K) 

TTTKFP. 
TOV/FR 
TPOL(I) 

VISFR(I,J) 

VISFRA 

VISFRE 

X(I),Y(I) 

XA(I),YA(I) 

XIC(I,J) 
YIC(I,J) 

XL,YL 


XLOC(I,J) 
YLOC(I,J) 


=  Time  to  fire  a  round  following  a  hit 
for  weapon  system  type  k  (k  =  1,2). 

=  Time  to  fire  first  round  after  target 
has  been  acquired  for  weapon  system 
type  k  (k  =  1,2). 

=  Time  to  fire  a  round  following  a  miss 
for  weapon  system  type  k  (k  =  1,2). 

=  Firing  rate  for  attacking  weapon  system, 

=  Firing  rate  for  defending  weapon  system. 

=  Total  percentage  of  lost  since  battle 
began  for  unit  i. 

=  The  fraction  of  unit  1  seen  by  unit  j. 

=  Fraction  of  unit  A  as  seen  by  unit  B. 

=  Fraction  of  unit  B  as  seen  by  unit  A. 

=  Coordinates  of  unit  i. 

=  Coordinates  of  alternate  position  for 
defender  i. 

=  Coordinates  of  the  j  ^  interval  endpoint 
of  the  route  for  unit  i. 

=  Distance  added  to  previous  interval 
endpoint  for  vehicle  to  move  DST  during 
a  time  step. 

th 
=  Coordinates  of  the  j    node  inputed  by 

the  user  for  the  route  of  unit  i. 
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